AD-A0B1  149 

[  UNCLASSIFIED 

ROYAL  AIRCRAFT  ESTABLISHMENT  FARN60R0U8H  (EN8LAND) 
COMPUTER  PROBRAMS  FOR  THE  DESI8N  OF  ACTIVE  RC  FILTERS 
MAR  79  LJ  RICHARDS 

RAE-TR-7903B  DRIC  -BR-70157 

,(U) 

F/8  9/S 

NL 

r 

lor| 

AO 

A08M49 

■ 

■ 

H 

H 

111113B1 

Hi 

nd  end  jp 

1 

1 

■ 

TR  79038 


o> < 

J. 

t— I  ! 


23 

Q- 

Ol 

3 

vjJ 


UNLIMITED 


TR  79038 


'JS® 


Bjc 


ROYAL  AIRCRAFT  ESTABLISHMENT 


Technical  Report  79038^ 

March  1979 

COMPUTER  PROGRAMS 
FOR  THE  DESIGN  OF 
ACTIVE  RC  FILTERS 


DTIO 

/r>  LECTEI 


Ft. 3  2  81980| 


D 


L.J.  Richards 


Procurement  Executive,  Ministry  of  Defence 
Farnborough,  Hants 


i  rr*U- 

I  :r. 


80  2  25  107  j 

"  . . .  Jd 


REPORT  DOCUMENTATION  PAGE 

Overall  security  classification  of  this  page 


UNCLASSIFIED 


As  far  as  possible  this  page  should  contain  only  unclassified  information.  If  it  is  necessary  to  enter  classified  information,  the  box 
above  must  be  marked  to  indicate  the  classification,  e.g.  Restricted,  Confidential  or  Secret. 


1.  DRIC  Reference  2.  Originator’s  Reference  3.  Agency  4.  Report  Security  Classification/Marking 

(to  be  added  by  DRIC)  Reference 


RAE  TR  79038 


UNCLASSIFIED 


5.  DRIC  Code  for  Originator  6.  Originator  (Corporate  Author)  Name  and  Location 

7673000W  Royal  Aircraft  Establishment,  Farnborough,  Hants,  UK 


5a.  Sponsoring  Agency’s  Code  6a.  Sponsoring  Agency  (Contract  Authority)  Name  and  Location 


7.  Title  Computer  programs  for  the  design  of  active  RC  filters 


7a.  (For  Translations)  Title  in  Foreign  Language 


7b.  (For  Conference  Papers)  Title,  Place  and  Date  of  Conference 


9b.  Authors  3, 4  .... 


13.  Project 


8.  Author  1 .  Surname,  Initials  9a.  Author  2 
Richards,  L.J. 


11.  Contract  Number  12.  Period 

N/A  N/A 


15.  Distribution  statement 

(a)  Controlled  by  - 

(b)  Special  limitations  (if  any)  - 


16.  Descriptors  (Keywords)  (Descriptors  marked  *  are  selected  from  TEST) 

Active  filters.  RC  filters 


10.  Date  Pages  Refs. 
March  90  I  1 3 
1979  1  V  1  J 


14.  Other  Reference  Nos. 
Rad-Nav  76 


17.  Abstract 

v 

This  Report  describes  three  Algol  68-R  programs  which  may  be  used  as  aids 
to  the  design  of  active  RC  filters.  Two  programs  calculate  the  performance  of 
standard  all-pole  filter  types,  the  first  program  calculating  the  attenuation  and 
the  second  the  group  delay.  The  third  calculates  component  values  needed  to 
implement  the  filter  as  a  cascade  connection  of  first  or  second  order  sections. 
Sensitivity  and  gain  performance  may  also  be  calculated.  Algorithms  are  described 
which  ensure  that  where  possible  component  values  lie  within  specified  bounds  and 
take  standard  values.  Specifications  of  the  programs  are  given,  together  with 
examples  of  their  use. 


ijj ,;i  it  /: 


UDC  621.372.852. 1 


ROYAL  AIRCRAFT  ESTABLISHMENT 


t  7  Technical  ^epi?t#^9038  |  I  ) 

Received  for  printing' 29  March  1979  — 


)  Alx  H  ! 


^COMPUTER  PROGRAMS  FOR  THE_DESIGN  OF  ^ACTIVE  RC  FILTERS < 

by 

j  L.  J. ^Richards  .  > 


hh 


jir\  f  fa  £  -TK-IV  i  y 


SUMMARY 


• — -^This  Report  describes  three  Algol  68-R  programs  which  may  be  used  as  aids 
to  the  design  of  active  RC  filters.  Two  programs  calculate  the  performance  of 
standard  all-pole  filter  types,  the  first  program  calculating  the  attenuation  and 
the  second  the  group  delay.  The  third  calculates  component  values  needed  to 
implement  the  filter  as  a  cascade  connection  of  first  or  second  order  sections. 
Sensitivity  and  gain  performance  may  also  be  calculated.  Algorithms  are 
described  which  ensure  that  where  possible  component  values  lie  within  specified 
bounds  and  take  standard  values.  Specifications  of  the  programs  are  given, 
together  with  examples  of  their  use.  I  -  — . . . 


Departmental  Reference:  Rad-Nav  76 


acaior 

~  J  iwuldt  1 
*1  -  .  _ 

-  -f.  A  ,:> 

* ■'  - c 1 1 -n 


Copyright 

© 

Controller  HMSO  London 
1979 


Thin  c  : 

for  p  j- 

diohibv  .'o:'  i  - 


s;t>o:  iifl 
h-r.itod. 


1  «••• : ;  a 

Lfl! 

j  i-i- eel 

’•-'''•/or 


i 


2 


LIST  OF  CONTENTS 

Page 

t  INTRODUCTION  3 

2  NORMALISATION  AND  DENORMALISATION  OF  TRANSFER  FUNCTIONS  4 

3  FILTER  TYPES  6 

3.1  Butterworth  filter  7 

3.2  Chebyshev  filter  8 

3.3  Gaussian  filter  8 

3.4  Bessel  filter  8 

3.5  Ultraspherical  family  of  filters  9 

3.6  Papoulis  filter  9 

3.7  Halpern  filter  9 

3.8  Synchronously  tuned  filter  10 

3.9  User-defined  filter  10 

4  CIRCUIT  TYPES  10 

4.1  First  order  section  11 

4.2  Unity  Gain  1 1 

4.3  Single  Feedback  12 

4.4  Multiple  Feedback  14 

4.5  State  Variable  16 

4.6  Ring  of  Three  18 

4.7  Mitra-Aatre  19 

4.8  Choice  of  circuit  types  20 

5  CHOICE  OF  COMPONENT  VALUES  21 

6  SENSITIVITY  23 

7  PERFORMANCE  OF  CIRCUIT  24 

8  GENERAL  NOTES  ON  THE  PROGRAMS  27 

9  PROGRAM  AC FA  30 

9. 1  Input  parameters  30 

9.2  Program  output  31 

9.3  Restrictions  31 

9.4  Examples  of  use  32 

10  PROGRAM  ACFD  32 

10. 1  Input  parameters  32 

10.2  Program  output  33 

10.3  Restrictions  33 

10.4  Example  of  use  34 

1 1  PROGRAM  ACFC  34 

11.1  Input  parameters  34 

11.2  Program  output  38 

11.3  Restrictions  40 

11.4  Examples  of  use  41 

1 2  CONCLUSIONS  42 

Appendix  A  Gain  and  frequency  transformations  43 

Appendix  B  Algorithm  for  satisfactory  component  values  45 

Appendix  C  Differential  sensitivity  47 

Appendix  D  Derivation  of  transmission  matrix  49 

List  of  symbols  and  abbreviations  53 

Illustrations  Figures  1-22 

Report  documentation  page  inside  back  cover 


038 


3 


n 

1  INTRODUCTION 

The  process  of  designing  an  analogue  filter  can  be  split  into  the  follow¬ 
ing  steps: 

j  (i)  Decide  on  the  mathematical  form  of  the  filter,  ie  compute  its 

transfer  function. 

(ii)  Decide  on  a  particular  type  of  circuit  to  realise  the  transfer 
function. 

(iii)  Compute  the  component  values  for  the  circuit. 

(iv)  Investigate  the  performance  of  the  circuit  thus  designed.  If  not 
satisfactory  repeat  one  or  more  of  the  previous  steps. 

None  of  these  steps  can  be  made  completely  automatic  but  the  use  of  computer 
programs  can  ease  the  designer's  task.  This  Report  describes  three  such  programs 
written  for  use  on  the  ICL  1906S  computer. 

Two  programs  calculate  attenuation  and  delay  performance  of  low-pass, 
high-pass  and  band-pass  filters  of  standard  all-pole  types  including  Butterworth, 
Chebyshev  and  Bessel  types.  The  first  program,  ACFA,  calculates  the  attenuation 
performance  and  the  second,  ACFD,  calculates  the  delay  performance.  The  third 
program,  ACFC,  allows  the  user  to  specify  one  of  a  number  of  standard  active  RC 
circuit  configurations  and  calculates  the  values  of  the  resistance  and  capacit¬ 
ance  values  required.  For  some  circuit  types  this  program  also  evaluates  the 
sensitivity  of  each  passive  component,  the  input  and  output  impedances  of  the 
circuit  and  the  filter  performance  with  non-ideal  active  elements. 

We  now  discuss  in  more  detail  the  filter  types  and  circuit  configurations 
chosen.  The  transfer  functions  considered  are  of  the  all-pole  type,  that  is  their 
transfer  functions  can  be  written  in  the  form 


H(s) 


Ks 


n  n-1 

s  +  an_]S  +  ...  -  a,s  +  aQ 


where  K,  aQ,  a(. 


a  .  are  real  constants  and  1=0  or  n  or  n/2  for 
n- 1 


low-pass  or  high-pass  or  band-pass  filters,  respectively.  Bandstop  and  all-pass 
filters  are  not  considered.  This  form  of  transfer  function  encompasses  for 
example  Butterworth  and  Chebyshev  filters  but  does  not  include  filters  such  as 
Elliptic  or  Inverse  Chebyshev  which  have  points  of  infinite  attenuation  at  finite 
frequencies . 
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The  filter  types  considered  are  listed  and  defined  in  section  3  below. 

The  common  practice  is  followed  of  defining  the  transfer  function  in  normalised 
form.  In  section  2  therefore  we  describe  the  normalisation  and  denormalisation 
procedures  and  discuss  the  frequency  transformations  used  to  obtain,  for  example, 
a  high-pass  filter  transfer  function  from  a  low-pass  prototype.  Also  in  section 
2  we  define  some  of  the  filter  parameters  used  in  the  rest  of  the  Report. 

The  circuits  considered  are  active  RC  circuits  using  operational  amplifiers 
as  the  active  elements.  Each  circuit  consists  of  a  number  of  subjircuits  or 
sections  connected  in  cascade.  Each  subcircuit  realises  a  first  or  second  order 
transfer  function.  Only  one  type  of  first  order  section  is  considered,  but  there 
is  a  choice  of  six  configurations  for  second  order  sections.  These  types  are 
listed  and  compared  in  section  4.  Although  different  circuits  may  use  different 
section  types  it  is  assumed  that  any  one  circuit  is  composed  of  only  one  section 
type. 

In  section  5  we  discuss  the  choice  of  component  values  and  explain  how  the 
computer  program  uses  preferred  capacitor  values  and  ensures  that  both  capacitor 
and  resistor  values  are  within  specified  bounds.  The  calculation  of  differential 
sensitivities  and  large  change  sensitivities  is  explained  in  section  6  and  in 
section  7  we  discuss  the  calculation  of  input  and  output  impedances,  voltage  range 
and  gain  performance  with  non-ideal  operational  amplifiers.  The  next  four  sec¬ 
tions  (8,  9,  10  and  11)  give  detailed  specifications  of  the  three  computer 
programs  and  give  examples  of  their  use. 

Active  filters  provide  an  alternative  to  traditional  passive  filters  for 
frequencies  up  to  a  few  hundred  kilohertz.  For  frequencies  up  to  a  few  kilohertz 
the  active  filter  becomes  particularly  attractive  since  it  is  more  amenable  to 
miniaturisation  than  its  passive  counterpart.  Active  filters  have  been  studied 

by  May  and  Brown1',  who  designed  a  low  frequency  narrow-band  filter  for  airborne 

12  ...  13 

use.  Kimbell  used  active  circuits  to  realise  pulse  forming  filters.  Norman 

considered  the  application  of  active  transitional  Butterworth-Thomson  filters  as 

pre-sampling  filters  in  telemetry  systems.  It  is  anticipated  that  active  filters 

will  continue  to  be  of  interest  to  RAE  for  low-frequency  applications,  for  example 

in  the  processing  of  voice  and  data  signals,  and  it  is  to  meet  this  requirement 

that  the  present  programme  of  work  has  been  undertaken. 

2  NORMALISATION  AND  DENORMALISATION  OF  TRANSFER  FUNCTIONS 

The  transfer  function  is  defined  to  be  the  ratio  of  the  Laplace  Transform 
of  the  output  to  that  of  the  input  and  we  denote  it  by  H(s)  .  It  is  convenient 
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to  write  an  nth  order  transfer  function  as  the  product  of  a  number  of  first  or 
second  order  transfer  functions 


H(s) 


Cks  +  Dks  +  Ek 
2 

s  +  A  s  +  B, 
k  k 


if  n  is  even,  or 


H(s)  = 


(n-1)/2  2 

ds  +  e  TT  CkS  +  °kS  +  Ek 


s  +  a 


n 

k=  1 


+  A  s  +  B, 
k  k 


if  n  is  odd. 

For  low-pass  filters  C^,  =  1^  =  d  =  0  (k  =  1,  2,  ...) 

For  high-pass  filters  Dk  =  Ek  =  e  =  0  . 

For  band-pass  filters  Ck=Ek=d=e=0  (only  even  orders  are  allowed). 

The  standard  filter  types  of  section  3  are  defined  only  for  low-pass 
'prototype'  filters  with  cut-off  frequency  1  rad/s  and  unity  gain.  We  can 
proceed  from  the  transfer  function  of  such  a  normalised  filter  to  obtain  the 
transfer  function  of  a  low-pass,  high-pass  or  band-pass  filter  of  arbitrary  cut¬ 
off  frequency  by  the  rules  given  in  Appendix  A. 

The  second  order  transfer  function  is  often  written  as 


Gw, 


0 


2  2 
s  +  UqS/Q  +  Wq 


for  a  low-pass  filter 


or 


Gs 


2  2 
s  +  WgS/Q  +  w 


or 


GWqS/Q 


2  2 
s  +  u>0s/Q  +  <oQ 


for  a  high-pass  filter 


for  a  band-pass  filter. 


where  G  is  the  gain.  The  quality  factor  Q  and  the  resonant  frequency  Wq 
are  related  to  the  A  and  B  coefficients  by 
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u0  =  ^ 

Q  =  /B/A  . 

Some  of  the  section  types  of  section  4  invert  the  output  with  respect  to 
the  input,  so  that  for  these  sections  G  should  strictly  be  negative.  In  this 
Report,  however,  G  is  assumed  to  be  positive,  ie  G  is  the  absolute  value  of 
the  gain.  Thus,  the  most  general  form  of  the  low-pass  transfer  function,  for 
example ,  is 


2  2 
s  +  WqS/Q  +  Wq 


Section  4.8  indicates  which  of  the  sections  are  inverting  and  which  non-inverting. 

Finally,  the  definition  of  the  cut-off  frequency  must  be  given.  For  filters 
having  a  monotonic  amplitude  response  in  the  pass-band  (eg  Butterworth,  Papoulis, 
Halpern,  Synchronously  Tuned)  it  is  conventional  that  the  cut-off  frequency  is 
that  frequency  at  which  the  amplitude  is  I//2  times  its  maximum  value,  ie  at 
which  the  attenuation  is  3  dB.  For  filters  having  an  oscillatory  amplitude  res¬ 
ponse  in  the  pass-band  (eg  Chebyshev,  Ultraspherical)  the  cut-off  frequency  is 

taken  to  be  the  highest  frequency  at  which  the  amplitude  differs  from  ics  maximum 
value  by  the  amount  of  the  ripple.  This  is  illustrated  in  Fig  8. 

3  FILTER  TYPES 

In  this  section  we  define  those  filter  types  used  in  the  programs.  The 
filter  types  are:  Butterworth,  Chebyshev,  Gaussian,  Bessel,  Ultraspherical, 
Papoulis,  Halpern  and  Synchronously  Tuned.  The  first  four  of  these  are  well- 
known^  and  details  of  all  the  types  are  given  in  Ref  2,  so  here  we  give  the  bare 
definitions  only.  In  each  case  we  give  the  definition  for  an  nth  order  low-pass 

filter  having  cut-off  frequency  1  rad/s  and  minimum  attenuation  0  dB.  For  typo¬ 

graphical  convenience  we  define  the  inverse  of  the  transfer  function,  ]/H(s)  , 

rather  than  the  transfer  function  H(s)  itself.  For  some  filters  it  is  more 

2 

convenient  to  quote  J H (s ) |  rather  than  H(s) .  To  obtain  the  transfer  function 
in  these  cases  one  solves  |h(s)|^  =  0  and  rejects  those  roots  with  positive 
real  parts.  The  remaining  roots  are  those  of  H(s)  =  0  . 

Once  one  has  the  transfer  function  it  is  straightforward  to  calculate  the 
gain  and  attenuation  at  any  frequency,  w  .  The  gain  is  simply 
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|H(s)  | 


where  s  =  jw  ,  and  the  attenuation  is 

-20  logl0|H(s)|  . 

The  group  delay,  D(io)  ,  defined  by 

0O)  =  -  arg  H  (s) 

is  slightly  more  complicated.  It  can  be  found  from  the  formula 


D(u)) 


n/2 

E 


k^1  <!k  -  “2>2  ♦  <v>2 


if  n  is  even,  or 


a  .  <nr/2  V\  +  -2) 


77'  E 


^  (Bk  -  W2)2  +  (A^)2 


if  n  is  odd. 

It  should  be  mentioned  here  that  not  all  the  filter  types  listed  above  are 
available  in  all  three  programs.  Papoulis,  Halpern  and  Synchronously  Tuned 
filters  are  not  available  in  either  of  programs  ACFC  (for  calculation  of  compon¬ 
ent  values)  or  ACFD  (for  calculation  of  delay).  This  is  not  too  great  a  restric¬ 
tion  on  the  usefulness  of  the  programs,  since  these  three  filter  types  are  rela¬ 
tively  rarely  used. 

3.1  Butterworth  filter 


n/2 

l/H(s)  =  1  [[s2  +  2  cos  ((2k  -  l)W(2n))s  +  ll 

k=i 


if  n  is  even,  or 


(n-l)/2 


TT  [• 


(s  +  1)  j  |  |s  +  2  cos(kir/n)s  +  1 
k=  I 


if  n  is  odd. 
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This  gives  an  attenuation  of  3  dB  at  w  =  1  . 
3 . 2  Chebyshev  filter 


n/2 

1/H(s)  =  ]j  (s  -  pk  +  jqk)(s  -  pk  -  jqk) 

k=  1 


if  n 

is  even,  or 

( s  +  sinh  (—  sinh  ^ 

\  Vn 

\ (n-l)/2 

0  17 

'  k=i 

jqkKs  "Pu¬ 

jqk> 

if  n 

is  odd. 

where 

.  ir  / 1  +  2k\ 
pk  =  '  Sln  2  {  n  ) 

sinh/i-  sinh  ' 

\n  t) 

TT  (l  +  2k\ 

%  '  “s  2  (,  „  ) 

cosh  (—  sinh  *  —  J 

Vn  e/ 

e2  -  to0'"1  -  I 

and  y 

is  the  attenuation  at 

to  =  1  . 

3.3 

Gaussian  filter  '■ 

1 / |H(s) 1 2  =  1 

0  2  4  3  6 

2  .  p  to  ,  p  to 

P“  2!  3! 

tt  2n 

+  +  p  u 
**•  n! 

> 

where 

p  is  a  constant  chosen 

so  that  the  attenuation  at  to  =  1 

is  3 

3.4 

Bessel  filter 

1/H(s)  =  Bn(D0s)/Bn(0) 


where  Bn  (x)  is  the  nth  order  Bessel  polynomial  which  can  be  defined  recursively 

by 

B0(x)  =  I 

Bj  (x)  =  x  +  I 

Bn(x)  =  (2n  -  l)Bn_j(x)  +  x2B^_2(x)  (n  5.  2)  , 


038 


pro 


9 


and  Dq  is  a  constant  chosen  so  that  Lhe  attenuation  at  is 


=  !  is  i  d B . 


3 . 5  Ultraspher ical  family  of  filters 


i  / !  h  ( s ) 


+  e 


Y  (x) 
n 


where  Y^(x)  is  a  polynomial  defined  recursively  by 


v*>  = 


Yj  (x)  =  x 


(2a  +  n)Y  (x) 
n 


and 


x(2a  +  2n 
2 


1  )  Y  (x) 

n- 1 


(n 


DY“  ,(x) 

n-z 


=  10°-IY  -  1 


(n  ^  2) 


where  y  is  the  attenuation  at  u  =  I  .  The  parameter  a  can  lie  in  the  range 
-1  <  a  <  +  00  .  Any  fixed  value  of  a  defines  one  filter  type  out  of  the  Ultra- 
spherical  family  of  filter  types. 


3.6  Papoulis  filter 


l/lH(s)j  =  1  +  T  (u) 


Tj  (to) 

= 

2 

u 

T2(u) 

= 

4 

w 

„  6 

4  2 

T3(oj) 

3w  - 

3io  +  io 

v 

3 
' — ✓ 

<3- 

H 

= 

6o8  - 

8iu8  +  3(o4 

(cu j 

= 

20w10 

-  40w8  + 

12 

10 

V“> 

50ui 

-  120u 

t7(«) 

= 

1 75u 1 4 

„r  12 
-  525a) 

V“> 

= 

490w 1 6 

-  1680w* 

6  o  4  *  2 

o  -  8w  +a) 

105<o8  -  40u^  +  6w4 
...  10  _  8  .  . 


)12  -  )624w’°  +  6 1 5oj8  -  1 20u;6  +  10u.4 


3. 7  Halpern  filter 


1  /  |  H(s)  i  2  =  1  +  Un(ui)  , 
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where  Uj (m) 

= 

2 

0)“ 

U,  (ul) 

= 

4 

L*} 

U3(u) 

= 

4mb  - 

4  2 

6u  +  3m 

U4(u>) 

= 

VO 

e 

CO 

1 

,  6  4 

1  6oj  +  8cj 

U5(u) 

= 

36m 1  ° 

-  90m8  +  80w6  -  30m4 

c  2 
+  5m 

U6  (w) 

= 

1  OOu) ' 2 

1 0  8 

-  288m  +  306m  -  1 

6  4 

144m  +  27m 

Uy  (ai) 

= 

14 

400oj 

-  I4u0m12  +  1932m10 

-  1330m8  +  476m' 

3.8  Synchronously  tuned  filter 

1 / | H(s) |2  =  [l  +  qw2] 

where  q  =  2  1  -  1 


3.9  User-defined  filter 

In  program  ACFC  (but  not  in  programs  ACFA,  ACFD)  the  user  may  define  his 
own  filter  type.  The  filter  is  specified  by  giving  the  poles  of  the  normalised 
low-pass  prototype  transfer  function.  It  is  assumed  that  all  the  poles  occur  in 
complex  conjugate  pairs,  with  the  possible  exception  of  a  single  real  pole. 

4  CIRCUIT  TYPES 

In  this  section  we  describe  the  conf igurat ions  used  to  realise  the  first 
and  second  order  transfer  functions.  For  each  circuit  type  we  give  the  design 
equations,  that  is,  equations  relating  the  component  values  to  terms  in  the 
transfer  function.  We  recall  that  a  first  order  transfer  function  has  the  form 


Ga 

s  +  a 

for  low-pass,  or 

Gs 

s  +  a 


for  high-pass  filter. 


A  second  order  transfer  function  has  the  form 


+ 


2 

s 


K _ 

+  As  +  B 
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for  a  low-pass,  or 


2 

s 


+ 


2 

US 

2 

s  +  As  +  B 


for  a  high-pass,  or 


2 

s 


+ 


Gu)qS/Q 


Q 


s  +  to 


2 

0 


2 

s 


Us 

+  As  +  B 


for  a  band-pass  filter. 

There  is  just  one  type  of  first  order  section  but  there  are  six  second 
order  sections:  Unity  Gain,  Single  Feedback,  Multiple  Feedback,  State  Variable, 
Ring  of  Three,  Mitra-Aatre.  These  circuits  are  just  a  few  out  of  the  many  that 
have  been  suggested  in  the  literature. 

In  each  case  the  transfer  function  is  realised  as  a  voltage  ratio.  On  the 
appropriate  circuit  diagram  the  input  voltage  is  denoted  by  V^n  and  the  output 
voltage  is  denoted  by  Vout  .  Where  the  same  configuration  can  be  used  for 
different  types  of  filter  VLp  denotes  the  low-pass  output,  VHp  denotes  the 
high-pass  output  and  Vgp  denotes  the  band-pass  output . 

3 

4 . 1  First  order  section 

This  is  illustrated  in  Fig  9  and  is  available  in  low-pass  and  high-pass 
forms.  The  design  equation  in  both  cases  is 


a  =  . 

Thus  C|  may  be  chosen  arbitrarily,  then  Rj  is  given  by 

Rj  =  I/CCja)  . 

The  gain  of  this  section  is  unity,  ie  G  =  1  . 

.  3 

4 . 2  Unity  Gain 

This  section  type  is  available  in  low-pass  or  high-pass  forms  and  is 
illustrated  in  Figs  10  and  11.  The  low-pass  configuration  uses  two  resistors, 
two  capacitors  and  two  unity  gain  operational  amplifiers.  The  design  equations 


are 
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A  =  1  /  (R ,  C  j  ) 

B  =  l/(RiR2C]C2) 
G  =  1  . 


and  C?  can  be  chosen  arbitrarily  and  Rj  and  R2  can  be  calculated  from 


Rj  =  1 / (AC  j ) 

R2  =  A/(BC2) 


The  high-pass  configuration  uses  two  capacitors,  two  resistors  and  one  unity 
gain  operational  amplifier.  The  two  capacitors  have  the  same  value.  The  design 
equations  are 


A  =  1 / (R2c2>  +  1 / ( R2  C , ) 

B  =  I /<R1R2CiC2) 

G  =  1  . 


We  can  choose  Cj  =  C2  arbitrarily,  then  put 

R,  =  A/(2RCj) 

R2  =  2/ (BCj ) 

4 

4 . 3  Single  Feedback 


This  circuit  type  is  available  in  low-pass,  high-pass  or  band-pass  forms 
as  shown  in  Figs  12,  13,  14. 


The  low-pass  configuration  has  four  resistors,  two  capacitors  and  one 
operational  amplifier.  The  design  equations  are 


B  =  l/(R3R4C2C3) 

g  -  Rjyy'Wi’ 


and  we  must  also  have 
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r3r4c3/(r]c1)  =  (R3  +  R4)R2/(Rj  +  r2)  . 

We  can  choose  C 2  arbitrarily,  then  put 

Rj  =  R2  =  x/(2Gy) 

R3  =  I / (yz) 

R4  =  (x  -  l/z)/y 
C  =  4GC2/x2 

C3  =  z2C2/(xz  -  1) 

where  y  =  C,^ 

x  =  A//B 


and  z  is  an  arbitrary  constant  satisfying  xz  >  1  . 

The  high-pass  configuration  has  three  resistors,  four  capacitors  and  one 
operational  amplifier.  The  design  equations  are 


A  = 


C4  (4  *  4) 


B  = 


and  we  must  also  have 


C  -  1-  R3)/(P2KjC3C4) 


Rl«2  *  R3)(CI  +  C2>  '  R2R3C4 


We  can  choose  C3  arbitrarily,  then  put 


Cj  =  c2  =  2gc3 


C4  =  z  C3/  (xz  -  I ) 


Rj  =  1 / (4Gxy) 

r2  =  i/(yz) 


R3  =  (x  -  l/z)/y 
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where  y  =  C3>/b 

x  =  A//B 

and  z  is  an  arbitrary  constant  satisfying  xz  >  1  . 

The  band-pass  configuration  uses  four  resistors,  four  capacitors  and  one 
amplifier.  Exact  design  equations  cannot  be  given  here  but  an  approximate  design 
procedure  is  to  choose  C2  arbitrarily  and  then  set 

ci  =  Gy/,“o 

C3  =  yz/((z  -  Do)0) 

C  =  —  z  r 

4  1  +  x  3 

Rj  =  1 / (Gy) 

R2  =  l/(yz) 

R3  =  (z  -  l)/(yz) 

R^  =  (x  +  DR^/z 

where  x  =  A//b 


y  =  c2m0/z 

and  z  is  an  arbitrary  parameter  such  that  z  >  1  . 

4 

4.4  Multiple  Feedback 

This  circuit  type  is  available  in  low-pass,  high-pass  and  band-pass  forms, 
as  shown  in  Figs  15,  16,  17. 

The  low-pass  configuration  uses  three  resistors,  two  capacitors  and  one 
operational  amplifier.  The  design  equations  are 


A  =  (1/Rj  +  1/R2  +  1/R3)/Cj 

B  =  l/(ClC2R2R3) 
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We  can  choose  arbitrarily  and  set 


Cj  =  4BC2(G  +  1)/A 
R3  =  A/(2BC2) 

R,  ■  R3/G 
r2  ■  R3/(G  ♦  1)  . 

The  high-pass  configuration  uses  two  resistors,  three  capacitors  and  one 
operational  amplifier.  The  design  equations  are 

A  =  (C,  +  C2  +  C3)/(R2C2C3) 

B  =  l/CRjR^C^) 

G  =  Cj/C3  . 


We  can  choose  Cj  =  C2  arbitrarily,  then  set 


C3  =  C,/G 

R]  =  A/(bC,[2  ♦  1/G]) 

R2  =  (2G  +  1 ) / (AC ) )  . 

The  band-pass  configuration  uses  three  resistors,  two  capacitors  and  one 
operational  amplifier.  The  design  equations  are 

A  =  (Cj  +  C2)/(R3CiC2) 

B  =  (1/Rj  +  1/R2)/R3 

G  =  Cj  /  (Rj  B)  . 

We  can  choose  Cj  =  C2  arbitrarily,  then  set 


Rj  =  1 / (CjGA) 


R2  =  x/(2/y  -  GA//B) 


R3  =  2x/y 
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where  x  =  1  /  (C  j  /b) 

y  =  A//B  . 

4.5  State  Variable^ 

The  single  circuit  shown  in  Fig  18  realises  a  low-pass,  high-pass  or  band¬ 
pass  transfer  function  depending  on  whether  V  ,  V  or  V  is  taken  as  the 

Lr  Hr  Hr 

output  voltage.  The  circuit  uses  six  resistors,  two  capacitors  and  three  opera¬ 
tional  amplifiers.  The  low-pass  design  equations  are 

At  =  (1  +  r2/r3>/ (i  +  VV 

Bt2  =  R2/R3 

G  =  (1  +  R3/R2)/(l  +  R4/R5) 

where  t  = 

We  can  choose  R3>  R^  and  Cj  arbitrarily.  Then,  if  we  let  x  =  R5/R4  * 
y  =  R2/R3  we  can  solve  the  following  quadratic  equation  for  Jy 

Q(G  -  l)y  +  /y  -  Q  =  0 


and  then  set 


x  =  QG>^y 

r2  =  R3y 

Rc  =  R. x 
5  4 

Rj  =  /y/CojgCj)  . 

However,  we  may  find  that  the  quadratic  equation  has  no  real  roots.  In  this  case 
we  abandon  the  attempt  to  achieve  the  specified  gain,  G  ,  and  set 

R2  =  R3 

R5  =  (2Q  -  1)R4 

R,  =  1/CwqCj)  . 

The  actual  gain  will  then  be  equal  to  2/(1  +  R^/R,.)  . 
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The  high-pass  design  equations  are 

At  «  (1  +  R2/R3)/(l  +  K5/R4) 

Bt2  =  R2/R3 

G  =  (!  +  R2/R3)/(l  +  R4/R5>  • 

We  can  choose  R3>  and  Cj  arbitrarily.  Then,  letting  x  =  R^/R^  » 

y  =  R2^R3  we  can  so^-ve  the  following  equation  for  /y 

Qy-^y+QCG+1)  =  0  . 

We  then  set 

x  =  QG/Sy 

R2  =  V 

Rc  =  R.  x 

R,  =  /  <co0c  j  )  . 

It  may  be  that  the  quadratic  equation  has  no  real  roots.  In  this  case  we  abandon 
the  attempt  to  achieve  the  specified  gain,  G  ,  and  set 

R2  =  R3 

R5  =  (2Q  -  1)R4 

ri  =  ‘/(“oV  ’ 

The  actual  gain  will  then  be  equal  to  2/(1  +  R^/R^). 

The  band-pass  design  equations  are 

At  =  (1  +  R2/R3)/(l  +  R5/r4) 

Bt2  =  R2/R3 

G  =  R5/R4  . 

We  can  choose  R3>  R^  and  Cj  arbitrarily.  Then  letting  x  =  R^/R^  * 
y  =  R2^R3  we  can  so^ve  c^e  following  quadratic  equation  for  /y 
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We  then  set 


Qy-(l+G)/y+Q  =  0 


x  =  G 


R2  ”  R3y 


R,  =  /y/CoipCj) 


It  may  be  that  the  quadratic  equation  has  no  real  roots, 
the  attempt  to  achieve  the  specified  gain,  G  ,  and  set 


In  this  case  we  abandon 


R5  =  (2Q  -  1)R4 

R,  =  l/O^Cj)  . 

The  actual  gain  will  then  be  2Q/(1  +  R^/R^ 

4 . 6  Ring  of  Three6  8 

The  circuit  of  Fig  19  realises  a  low-pass  or  band-pass  filter  depending  on 

whether  V  or  V  is  taken  as  the  output.  It  uses  six  resistors,  two 
LP  BP 

capacitors  and  three  amplifiers.  The  low-pass  design  equations  are 

“0  ■  V'VjVi'V 

Q  .  R5CA 

G  -  R4/R6  . 

We  can  choose  Rj ,  and  Cj  =  arbitrarily,  then  set 


Clu0 v  Rl 


R5  =  Q/(C]Wo) 
R6  =  R3/G 


R4  “  R3  • 
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The  band-pass  design  equations  are 


"0  "  R2/(R.R3R4C1C2) 


R5C2“o 


G  =  R5/R6 


We  can  choose  R,,  and  arbitrarily,  then  set 


R5  =  Q/CCjUq) 

R6  "  VG 

R4  =  R3  * 


,  9 

4.7  Mitra-Aatre 


This  circuit  is  available  only  in  band-pass  form  and  is  illustrated  in 
Fig  20.  It  uses  no  capacitors,  seven  resistors  and  two  operational  amplifiers. 
This  circuit  differs  from  all  others  considered  in  this  Report  in  that  it  uses  a 
different  model  of  the  operational  amplifier.  Other  circuits  assume  an  ideal 
amplifier,  ie  if  Vj (s)  and  V2(s)  are  the  inverting  and  non-inverting  inputs 
to  the  amplifier  at  a  given  frequency  and  V^Cs)  is  the  corresponding  output 
then 

V3(s)  =  Aq(v2(s)  -  V,  (s>) 


where  Aq  =  +  °°  . 

The  present  circuit  replaces  this  equation  by 


V3(s)  =  7  (V2(s)  "  Vs)) 

where  P  is  the  gain-bandwidth  product  of  the  amplifier.  Like  the  first  equa¬ 
tion  this  is  only  an  approximation  to  the  truth,  but  Mitra  and  Aatre  claim 
reasonable  performance  for  centre  frequencies  between  5  and  200  kHz. 

The  design  equations  are 
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B  =  P,  P2/  [(1  +  R3/R4)(l  +  r6/r7>] 
A  =  Pj/O  +  R5/Rj  +  R5/R2) 


We  can  choose  R^  and  R^  arbitrarily,  then  set 


r5  =  R!g 


where  Pj  and  P2  are  the  gain-bandwidth  products  of  the  two  amplifiers.  For 
ease  of  programming  we  have  assumed  that  all  operational  amplifiers  have  the  same 
gain-bandwidth  product,  but  this  is  not  essential. 


4 . 8  Choice  of  circuit  types 

Since  we  have  described  six  alternative  second  order  circuit  types  we  may 
ask  which  of  these  types  is  best.  No  definitive  answer  can  be  given  to  this 
question  since  each  type  has  certain  advantages.  The  Unity  Gain  circuit  has  the 
smallest  number  of  passive  components  and  has  a  simple  structure  but  cannot  give 
arbitrary  gain  and  is  not  available  in  band-pass  form.  Both  Single  Feedback  and 
Multiple  Feedback  circuits  use  only  one  amplifier  per  section  and  can  realise 
arbitrary  gain.  Of  the  two,  the  Multiple  Feedback  uses  fewer  passive  components, 
but  may  require  a  larger  spread  of  element  values  in  some  circumstances.  The 
Single  Feedback  circuit  may  be  viewed  with  some  suspicion  since  it  relies  on 
exact  pole-zero  cancellation  in  the  transfer  function.  Both  Single  Feedback  and 
Multiple  Feedback  types  use  in  general  at  least  one  capacitor  which  cannot  be 
constrained  to  take  a  standard  value. 

The  State  Variable  and  Ring  of  Three  circuits  both  use  three  amplifiers 
per  section.  They  use  the  same  number  of  passive  components.  The  Ring  of  Three 
circuit  can  achieve  arbitrary  gain,  whereas  the  State  Variable  circuit  can  only 
do  this  in  some  circumstances. 

The  Mitra-Aatre  circuit  uses  no  capacitors  but  requires  a  knowledge  of  the 
gain-bandwidth  product  of  the  amplifiers.  It  is  intended  for  use  at  higher 
frequencies,  say  greater  than  5  kHz. 
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If  a  rule  of  thumb  is  desired  we  may  say:  for  low  Q  sections  use  the 
Unity  Gain  circuit.  For  good  sensitivity  performance  use  the  Ring  of  Three 
circuit.  For  high  frequency-operation  try  the  Mitra-Aatre  circuit. 

Program  ACFC  assumes  that  all  second  order  sections  will  be  of  the  same 
type,  ie  either  all  Unity  Gain  or  all  Single  Feedback  etc.  But  there  is  nothing 
to  stop  the  designer  mixing  section  types  in  the  real  circuit.  For  example,  a 
sixth  order  filter  might  have  two  low  Q  sections  and  one  high  Q  section.  One 
could  run  program  ACFC  with  two  sets  of  data,  one  for  Unity  Gain  sections  only 
and  one  for  Ring  of  Three  sections  only.  One  could  then  use  the  Unity  Gain  com¬ 
ponents  for  the  lew  Q  sections  and  the  Ring  of  Three  components  for  the  high  Q 
section. 

As  we  mentioned  in  section  2,  some  section  types  invert  their  output  with 
respect  to  the  input.  The  inverting  types  are:  all  single  feedback,  multiple 
feedback.  Ring  of  Three,  Mitra-Aatre  types  and  the  high-pass  State  Variable. 

The  remaining  types  are  non-inverting,  ie  first  order  section,  Unity  Gain,  low- 
pass  and  high-pass  State  Variable. 

5  CHOICE  OF  COMPONENT  VALUES 

In  all  the  circuits  of  section  4  the  number  of  passive  elements  exceeds 
the  number  of  design  equations,  so  that  one  or  more  components  may  be  chosen 
arbitrarily.  We  can  use  this  freedom  in  order  to  achieve  satisfactory  values 
for  all  components.  We  would  like  the  component  values  to  have  the  following 
properties : 

(i)  All  capacitance  and  resistance  values  should  lie  between  lower  and 
upper  bounds  specified  by  the  user. 

(ii)  Only  preferred  capacitance  values  should  be  used. 

(iii)  Values  should  lie  near  the  middle  of  the  ranges  specified  by  the  user. 

(iv)  The  spread  of  resistance  values  and  the  spread  of  capacitance  values 
should  be  small. 

It  will  not,  in  general,  be  possible  to  achieve  all  of  these  properties: 
sometimes  we  shall  not  be  able  to  achieve  any  of  them.  Program  ACFC  tries  to 
obtain  satisfactory  component  values  by  an  iterative  process  that  we  now  describe. 

In  program  ACFC  arbitrary  capacitances  are  allowed  to  take  only  the  follow¬ 
ing  values  (in  farads) 
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10 

10 

10 

10 

10 


-11 

-10 

-9 

-8 

-7 


2.2  x 
2.2  x 

2.2  x 

2.2  x 


10 


■1  1 


10 


-10 


10 


-9 


10 


-8 


5  x  10 


-l : 


5  x  10 


-10 


5  x  10 


-9 


5  x  10 


These  preferred  values  could  of  course  be  changed  as  desired.  For  most  of  the 
circuits  (First  Order,  Unity  Gain,  State  Variable,  Ring  of  Three,  Mitra-Aatre) 
all  capacitances  are  arbitrary,  but  in  some  (single  feedback  and  multiple  feed¬ 
back)  there  are  non-arbitrary  capacitances  which  cannot  be  forced  to  have  pre¬ 
ferred  values. 

The  program  makes  no  attempt  to  use  only  preferred  resistance  values,  but 

it  does  allow  the  user  to  specify  lower  and  upper  bounds,  R  .  and  R  for 

*  '  min  max 

resistances,  and  also  C  .  and  C  for  capacitances.  The  user  can  also 

min  max  r 

provide  a  single  preferred  value  for  any  arbitrary  resistance,  though  if  this 
choice  leads  to  other  components  taking  unsatisfactory  values  it  may  be  changed. 

For  any  first  or  second  order  section  the  following  procedure  is  carried 

out: 

(i)  Set  all  arbitrary  capacitances  to  a  preferred  value  in  the  middle 

-9 

of  the  range,  namely  10  farads. 

(ii)  Set  all  arbitrary  resistances  to  the  preferred  value  supplied  by  the 
user. 

(iii)  Calculate  the  other  component  values  from  the  equations  given  in 
section  4. 

(iv)  Check  whether  any  calculated  component  values  lie  outside  their 
bounds.  If  so  change  one  of  the  arbitrary  values.  If  an  arbitrary 
capacitance  is  to  be  changed,  select  the  next  lower  or  next  higher 
preferred  value. 

(v)  Go  to  step  (iii). 

The  procedure  is  continued  until  either  a  satisfactory  set  of  components  is 
obtained  or  it  is  clear  that  no  such  set  can  be  obtained.  If  a  satisfactory  set 
cannot  be  found  an  appropriate  message  is  printed,  but  the  program  does  not  halt: 
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it  uses  the  last  calculated  set  of  values.  Having  calculated  the  component 
values  we  may  wish  to  change  them  prior  to  performing  a  sensitivity  or  performance 
analysis,  for  example  if  we  wish  to  investigate  the  effect  of  using  only  standard 
value  resistances.  The  program  allows  the  calculated  values  to  be  overwritten 
in  this  way. 

We  have  not  yet  described  in  detail  how  arbitrary  component  values  are 
changed  to  avoid  having  components  lying  outside  their  bounds.  Appendix  B  gives 
the  algorithm  used  for  the  Ring  of  Three  low-pass  circuit.  Similar  algorithms 
have  been  derived  for  the  other  circuit  types. 

6  SENSITIVITY 

It  is  of  interest  to  determine  how  the  performance  of  a  circuit  changes  as 
the  value  of  some  component  changes.  Let  r  denote  the  variable  measuring  the 
performance  and  p  denote  the  parameter  being  changed.  We  define  the  differen¬ 
tial  relative  sensitivity  to  be 

Sr  =  lim  (Ar/r) /(Ap/p) 
p  Ap-K) 

£ 

r  Sp 

We  are  mostly  interested  in  the  gain  of  the  filter  which  is  equal  to  the  absolute 
value  of  the  transfer  function.  It  is  shown  in  Appendix  C  that  the  required 
sensitivity  can  be  expressed  in  terms  of  the  quantities 

n  3U>0  9Q  9G 

P,VQ,G’  9P  •  9P  •  aF  • 

To  evaluate  the  partial  derivatives  we  need  to  consider  the  design  equations  of 
the  particular  circuit.  For  example,  for  the  Ring  of  Three  circuit, 


2 

W0  " 

'j'd’iWiS* 

30 

3u0 

9RI  = 

-  0.5u)q/Rj 

3“o 

3R2 

0. 5Uq/R2 

etc. 
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Having  found  the  differential  sensitivities  we  can  estimate  the  tolerances 
required  for  each  passive  component.  For  suppose  that  M  is  the  maximum  allow¬ 
able  relative  change  in  gain,  ic  we  require  that  ■$  M  .  Then  since 


W  ap 


we  have 


<  M  , 

p  '  ' 


:so  tint  MVS'  is  an  esr.  imate  for  the  required  tolerance  of  parameter  p  . 

I  H  [ 

The  differential  sensitivity  S  1  and  the  corresponding  tolerance  figure 

P 

give  a  useful  indication  of  how  the  performance  will  be  changed  by  small  changes 
in  the  parameter.  However  these  figures  should  not  be  given  too  much  importance 
since  they  assume  inf initesimal  changes  in  the  parameters  and  also  assume  that 
only  one  parameter  is  changing  at  a  time,  thus  ignoring  coupling  effects  between 
components.  For  these  reas%  as  it  is  worth  calculating  also  large  change  sensi¬ 
tivities,  by  which  we  mean  relative  changes  in  gain  produced  by  specified  changes 
in  one  or  more  components.  The  computer  program  allows  the  user  to  increase  or 
decrease  any  or  all  of  the  components  by  arbitrary  amounts  and  calculates  the 
corresponding  gain  from  the  design  equations  of  the  appropriate  section  type. 

The  program  does  not  however  perform  any  automatic  or  ’Monte  Carlo’  calculations 
to  establish  the  effect  of  random  perturbations  of  the  parameters.  It  is  left  to 
the  user  to  decide  which  sets  of  parameter  changes  to  investigate. 

Program  ACFC  performs  sensitivity  calculations  for  only  two  filter  types 
(Unity  Gain  and  Ring  of  Three) ,  but  there  is  no  reason  in  principle  why  other 
filter  types  should  not  be  included. 

7  PERFORMANCE  OF  CIRCUIT 

If  the  design  has  been  correctly  carried  out  it  would  seem  that  the  circuit 
should  behave  exactly  as  desired.  There  are  a  number  of  reasons  why  this  might 
not  be  true.  The  component  values  as  calculated  from  the  equations  of  section  4 
may  have  been  overwritten  as  described  in  section  5,  for  example  if  resistances 
have  been  rounded  to  convenient  values.  Also,  the  operational  amplifiers  may 
have  been  incorrectly  modelled.  The  equations  of  section  4  assume  that  the  opera¬ 
tional  amplifiers  have  infinite  input  impedances  and  zero  output  impedances.  The 
equations  for  all  section  types  except  the  Mitra-Aatre  assume  the  amplifiers  have 
infinite  open-loop  gain  at  all  frequencies.  Program  ACFC  allows  the  user  to 
specify  values  of  input  and  output  impedances  and  open-loop  gain,  the  same  values 
being  assumed  to  apply  to  all  the  amplifiers  in  the  circuit.  The  program  will 
then  calculate,  for  any  frequency  specified  by  the  user, 
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(i)  the  input  and  output  impedances  and  gain  of  each  section, 

(ii)  the  minimum  and  maximum  voltages  occurring  in  the  circuit,  and  the 
voltage  spread,  (ie  the  ratio  of  maximum  to  minimum  voltage), 

(iii)  the  transmission  matrix  of  the  complete  circuit, 

(iv)  the  input  and  output  impedances  and  gain  of  the  complete  circuit. 


This  should  help  the  user  to  make  a  good  estimate  of  the  circuit  performance,  but 
it  will  not  enable  him  to  predict  exactly  the  actual  circuit  performance  since 
there  are  many  factors  not  taken  into  account  by  the  program.  For  example,  even 
with  realistic  input  and  output  impedances  and  open-loop  gain  the  model  of  the 
operational  amplifier  is  crude  -  it  does  not  include  effects  due  to  frequency  - 
dependent  gain,  slew  rate,  output  clipping  etc. 


We  now  discuss  the  calculation  of  the  quantities  (i)  to  (iv)  above.  The 
zero  output  impedance  of  the  ideal  operational  amplifier  means  that  each  section 
can  be  considered  independently  of  the  rest  of  the  circuit.  When  non-ideal 
amplifiers  are  considered  however,  the  different  sections  will  interact.  One 
consequence  of  this  is  that  the  ordering  of  the  sections  makes  a  difference  and 
must  be  specified.  If  the  user  is  not  sure  which  ordering  to  use  he  may  either 
try  different  orderings  and  see  which  gives  the  best  results  or  use  the  rule  of 
thumb  that  the  sections  should  be  arranged  in  order  of  increasing  Q  .  Another 
consequence  is  that  we  are  unable  to  derive  a  voltage  transfer  function  since  we 
do  not  know  how  the  network  is  terminated.  The  best  we  can  do  is  to  derive  the 
transmission  matrix*^  of  the  network,  ie  the  matrix  relating  input  to  output 
quantities  as 
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where  suffix  1  denotes  input  quantities  and  suffix  2  denotes  output  quantities 
(the  positive  inwards  current  convention  is  assumed  here).  Each  section  has  an 
associated  transmission  matrix  and  since  the  sections  are  connected  in  cascade, 
the  transmission  matrix  of  the  network  is  the  product  of  the  transmission  matrices 
of  the  individual  sections.  Once  we  have  the  transmission  matrix  it  is  easy  to 
calculate  the  input  and  output  impedances  and  gain.  If 


1 1 

3 1 2 

'21 

322 

26 


then 


and 


the  input  impedance 


v° 


the  output  impedance 


=0 


the  gain 


V0 


(NB.  The  transmission  matrix  elements  are  usually  denoted  by  A,  B,  C,  D  where 
A  =  a  j,  B=aj2»  C  =  a^j,  D=  a22‘  ^ave  av°ided  this  notation  since  we  have 
reserved  the  letters  A,  B,  C,  D  for  other  uses). 

The  voltage  spread  can  be  estimated  as  follows.  The  voltage  and  current 
at  the  end  of  the  last  section  are  set  arbitrarily  to  1  and  0  respectively. 
Pre-multiplication  by  the  transmission  matrix  of  the  last  section  gives  tl.e 
voltage  and  current  at  the  beginning  of  the  last  section.  Multiplication  by  each 
matrix  in  turn  thus  gives  the  voltages  at  the  end  of  each  section.  We  can  then 
pick  out  the  largest  and  smallest  of  these  voltages  and  divide  the  largest  by  the 
smallest  to  give  the  voltage  spread.  It  should  be  remembered  that  these  figures 
are  based  on  unit  voltage  and  zero  current  at  the  end  of  the  last  section.  It 
should  also  be  noted  that  the  voltages  inside  each  section  are  not  calculated. 

It  is  possible  that  such  a  voltage  is  greater  than  the  largest  of  the  inter¬ 
section  voltages.  Thus  the  voltage  spread  as  calculated  is  only  a  lower  bound 
for  the  true  voltage  spread. 

The  problem  of  finding  the  transmission  matrix  of  a  single  section  is 
discussed  in  Appendix  D.  The  calculations  will  depend  on  the  section  type,  but 
as  an  example,  the  Appendix  gives  the  procedure  for  the  low-pass  Ring  of  Three 
section. 

So,  having  determined  the  circuit  component  values,  quantities  (i)  to  (iv) 
can  be  found.  This  assumes  that  we  have  decided  on  the  order  in  which  the  sec¬ 
tions  are  connected.  For  the  calculation  of  component  values  and  sensitivities 
the  ordering  of  the  sections  is  irrelevant  since  the  operational  amplifiers  are 
assumed  to  be  ideal  for  these  calculations.  When  realistic  input  and  output 
impedances  are  used,  however,  the  ordering  does  make  a  difference,  and  must  be 
specified.  There  is  no  simple  procedure  for  choosing  the  best  ordering,  but  a 


27 


rule  of  thumb  that  is  sometimes  used  is  to  arrange  the  sections  in  order  of 
increasing  Q  .  Program  ACFC  allows  the  user  to  specify  the  ordering  or  to  use 
this  rule  of  thumb. 

In  this  section,  we  have  described  how  various  performance  characteristics 
of  the  circuit  (items  (i)  to  (iv)  above)  can  be  calculated  for  any  circuit  type. 
Program  ACFC  allows  these  characteristics  to  be  calculated  for  only  two  circuit 
types  (Unity  Gain  and  Ring  of  Three),  but  there  is  no  reason  in  principle  why 
the  same  calculations  should  not  be  performed  for  any  other  circuit  type. 

8  GENERAL  NOTES  ON  THE  PROGRAMS 

In  this  section  we  give  notes  on  programs  ACFA,  ACFD  and  ACFC  and  on  their 
use  on  the  RAE  ICL  1906S  computer. 

The  programs  are  written  in  Algol  68-R  and  are  designed  to  be  run  under  the 
George  4  operating  system.  The  following  files  are  contained  in  directory 
:RN4660G. 

ACFA,  ACFC,  ACFD,  ACFSEG1 ,  ACFSEG2 ,  ACFSEG3,  ACFALB,  ACFALBUM,  ACFTIDY. 

Files  ACFA,  ACFC  and  ACFD  contain  source  programs.  ACFSEG1,  ACFSEG2,  ACFSEG3 
contain  source  procedures  used  by  these  programs.  Files  ACFALB  and  ACFALBUM  are 
Algol  68-R  albums  containing  compiled  versions  of  the  procedures  in  ACFSEG1 , 
ACFSEG2  and  ACFSEG3.  ACFTIDY  is  a  macro  used  to  edit  the  output  produced  by 
ACFC.  Copies  of  card  listings  of  these  files  may  be  obtained  from  the  author. 

The  three  programs  may  be  run  in  either  interactive  or  non-interactive  mode. 
In  interactive  mode,  the  programs  are  run  from  a  MOP  console  and  the  user  inputs 
data  in  response  to  questions  output  by  the  programs  to  the  console.  In  non¬ 
interactive  mode  the  data  are  first  stored  in  a  file,  the  program  is  then  run 
either  at  a  MOP  console  or  as  a  background  job,  the  results  being  stored  in  a 
file. 

The  commands  for  running  program  ACFA  in  non-interactive  mode  are 

COMPILE  : RN4660G.  ACFA 
AS  *TRO,  file  1 
AS  *LPO,  file  2 
EXECUTE 

where  file  I  is  the  file  containing  the  data  and  file  2  is  the  file  which  is  to 
contain  the  output.  If  running  in  interactive  mode  the  two  assignment  commands 
are  omitted,  so  the  commands  are  simply 
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COMPILE  :RN4660G.  ACFA 
EXECUTE. 

The  commands  for  running  ACFD  are  similar.  They  are 


COMPILE  :RN4660G.  ACFD 
AS  *TRO,  file  1 
AS  *LPO,  file  2 
EXECUTE . 


(if  in  non-interactive  mode) 


The  commands  for  running  ACFC  are  similar  except  that  an  extra  command  is  needed 
after  compilation 


COMPILE  :RN4660G.  ACFC 
GETEBM  PR0GC 
AS  *TRO,  file  T 
AS  *LPO,  file  2 
EXECUTE . 


(if  in  non-interactive  mode) 


If  program  ACFC  is  run  in  non-interactive  mode  the  output  file  will  contain  the 
questions  which  would  have  been  asked  at  the  console  if  the  mode  were  interactive. 
Since  these  questions  are  of  no  interest  to  the  user  after  the  program  has  been 
run,  a  macro,  ACFTIDY,  has  been  written  to  edit  out  the  questions  from  the  output 
file.  The  macro  command  is 


ACFTIDY  file  2 

if  the  job  is  being  run  under  directory  :RN4660G  or 

OBEY  :RN4660G.  ACFTIDY,  PARAM  (file  2) 

if  the  job  is  being  run  under  another  directory.  The  macro  will  generate  two 
files  named  LJRTEM1  and  LJRTEM2  and  then  erase  them. 

Normally  the  programs  will  read  all  their  data  and  calculate  all  the 
required  results.  They  may  fail  however  if  either  the  restrictions  given  in 
sections  9.3,  10.3,  11.3  are  not  observed,  or  if  the  transfer  function  cannot  be 
obtained.  In  a  number  of  cases  to  obtain  the  transfer  function  the  program  needs 
to  find  the  roots  of  a  polynomial  equation.  If,  particularly  for  higher  orders, 
these  roots  cannot  be  found  an  appropriate  message  is  printed  and  the  program 
will  be  halted. 


Program  ACFD  automatically  calls  up  and  makes  use  of  album  ACFALB  and  pro¬ 
gram  ACFC  uses  the  two  albums  ACFALB  and  ACFALBUM.  The  user  need  not  normally 
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concern  himself  with  these  two  albums  or  with  the  source  files  ACFSEG1 ,  ACFSEG2, 
ACFSEG3  from  which  they  are  obtained.  However,  for  reference  we  give  the  proce¬ 
dure  for  setting  up  the  two  albums  from  the  three  source  files.  To  set  up 
ACFALB  the  following  commands  are  required 


n 


CE  ACFALB  <*ED,  BUCKET!,  KW9RDS  19) 

COMPILE  ACFSEG1 
UPDATEALBUM  ACFALB 
CLEART6  :0 
PUTIN  SEG 1 
END 

COMPILE  ACFSEG2 
UPDATEALBUM  ACFALB 
PUTIN  SEG2 
END. 

To  set  up  ACFALBUM  one  should  first  set  up  ACFALB  and  then  issue  the  following 
commands 

CE  ACFALBUM  (*ED,  BUCKET  1,  KW0RDS  20) 

COMPILE  ACFSEG3 
UPDATEALBUM  ACFALBUM 
CLEART6  : ACFALB 
PUTIN  SEG3 
END. 

Also  for  reference,  we  point  out  that  if  the  albums  are  set  up  under  a 
directory  other  than  :RN4660G  then  the  first  line  of  some  of  the  files  will  need 
to  be  changed.  For  example  the  first  line  of  file  ACFD  is  at  present 

PR0G1B  'WITH'  SEG1 ,  SEG2  ’FROM’  :RN4660G.  ACFALB. 

If  the  albums  to  be  used  are  :ABC!23.  ACFALB  and  :ABC123.  ACFALBUM  then  the  first 
line  of  ACFD  should  be  changed  to 

PROG  1 B  'WITH'  SEG1 ,  SEG2  'FROM'  :ABC123.  ACFALB. 

Similar  changes  will  be  required  for  files  ACFC,  ACFSEG2  and  ACFSEG3. 

For  users  more  familiar  with  FORTRAN  than  Algol  68-R  it  should  be  observed 
that  in  input  and  output  the  symbol  denoting  exponent  of  10  is  the  ampersand  (&) 
rather  than  the  E  of  FORTRAN.  Thus  the  number  0.0012  could  be  written  as 
I.2E-3  in  FORTRAN  but  would  be  1.2&-3  in  Algol  68-R.  It  should  also  be  noted 
that  a  real  number  should  contain  a  decimal  point  and  at  least  one  digit 
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following  the  decimal  point.  In  FORTRAN  a  number  could  be  input  in  the  form 
1 2E-4  or  12  .E-R  but  in  A  '•  go  i  68-R  the  number  would  have  to  be  given  as  12.0&-4. 

On  output  it  should  h  r ...■  t.<:d  that  Algo]  68-R  outputs  a  complex  number  in  the  form 

real  part  ?  imaginary  part. 

9  PROGRAM  AC FA 

This  program  calculates  the  attenuation  performance  of  a  given  filter.  The 
filter  type  is  chosen  from  tile  list  given  in  section  3.  The  user  can  either 
*pecifv  the  order  of  the  filter  and  let  the  program  calculate  the  attenuation  at 
specified  frequencies  of  interest  or  let  the  program  calculate  the  minimum  order 
needed  to  achieve  a  specified  attenuation  at  a  specified  frequency. 

9. 1  Input  parameters 

The  following  input  parameters  are  required  in  the  following  order.  Unless 
stated  otherwise  there  should  be  one  item  per  record. 

(1)  One  of  the  three  symbols  L,  H,  B,  depending  on  whether  the  filter  is 
low-pass,  high-pass  or  band-pass.  The  symbol  should  be  the  first  character 
of  the  record. 

(2)  Filter  name.  This  should  begin  at  the  first  character  of  the  record. 
Only  the  first  four  characters  of  the  name  are  significant.  Thus  a  Gaussian 
filter  could  be  specified  by  either  GAUSSIAN  or  GAUS. 

(3)  Alpha.  This  real  number  applies  only  to  Ultraspherical  filters  (see 
section  3.5)  and  selects  one  filter  type  out  of  the  Ultraspherical  family 
of  filter  types.  For  other  filter  types  this  parameter  should  be  omitted. 

(4)  The  order  of  the  filter.  If  this  is  given  as  zero  the  minimum  order 
required  will  oe  calculated  by  the  program.  In  the  case  of  band-pass 
filters  the  or^er  to  be  given  is  the  order  of  the  low-pass  prototype 

(i<2  half  the  actual  order). 

(5)  The  attenuation  in  dB  required  at  a  given  stop-band  frequency.  This 
parameter  should  only  be  given  if  the  filter  order  was  given  as  zero. 

(6)  The  cut-off  frequency  in  hertz.  This  parameter  should  only  be  given 
if  the  filter  is  iow-pass  or  high-pass. 

(7)  Geometric  centre  frequency  and  bandwidth,  both  in  hertz.  These 
numbers  should  appear  on  the  same  record.  The  parameters  should  only  be 
given  if  the  filter  is  band-pass. 
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(8)  The  attenuation  in  dB  at  the  cut-off  frequency  (or  frequencies). 

(9)  The  number  of  frequencies  of  interest.  This  is  an  integer  not  greater 

than  20. 

(10)  The  frequencies  of  interest,  in  hertz.  If  the  order  of  the  filter 

was  given  as  zero  the  first  frequency  given  will  be  assumed  to  be  the 

frequency  at  which  the  attenuation  was  specified.  The  numbers  may  be 
spread  over  any  number  of  records. 

9.2  Program  output 

The  attenuation  at  each  frequency  of  interest  is  output  in  the  form 

ATTENUATION  AT  FREQ.  no.  WITH  ORDER  no.  IS  no. 

The  order  referred  to  is  the  order  of  the  low-pass  prototype  which  in  the  case 
of  band-pass  filters  is  half  the  actual  order.  To  emphasize  this,  in  the  band¬ 
pass  case  the  actual  order  of  the  band-pass  filter  is  output  after  the  last 
frequency  of  interest.  If  the  program  is  required  to  find  the  minimum  order 
needed  to  achieve  a  given  attenuation  at  a  given  frequency  then  the  attenuation 
at  all  frequencies  of  interest  will  be  output  for  orders  1,  2,  3  ...,  the  order 
being  increased  until  either  the  required  attenuation  is  achieved  or  the  order 
exceeds  the  maximum  order  allowable  for  the  filter  type  (see  section  9.3  below). 

9.3  Restrictions 

The  order  of  the  filter  types  are  restricted  as  shown  in  the  following 

table. 


Filter  type 

Butterworth 

Chebyshev 

Ultraspherical 

Papoulis 

Halpern 

Gaussian 

Bessel 

Synchronously  Tuned 


Maximum  order 
30 
20 
20 
8 
7 

20 

6 

20 


In  the  case  of  band-pass  filters  the  order  referred  to  in  this  table  is  that  of 
the  low-pass  prototype. 


As  stated  previously  there  can  be  not  more  than  20  frequencies  of  interest. 
If  it  is  desired  to  change  this  limit  to,  say,  50,  then  the  tenth  line  of  program 
ACFA  should  be  changed  from 


[l  :  20]  ’REAL'  WINT  ; 
to 

[1  :  50]  'REAL*  WINT  ; 

9.4  Examples  of  use 

(1)  A  fourth  order  high-pass  Bessel  filter  with  cut-off  frequency 
1500  Hz  is  to  be  considered.  The  attenuations  at  frequencies  1200  Hz, 

1000  Hz  and  800  Hz  are  required.  Fig  1  shows  the  results  of  an  inter¬ 
active  run  of  program  ACFA. 

(2)  It  is  required  to  construct  a  Chebyshev  band-pass  filter  with  centre 
frequency  500  Hz  and  bandwidth  5  Hz.  A  ripple  of  0.1  dB  is  allowed  in 
the  pass-band.  An  attenuation  of  at  least  70  dB  is  required  at  400  Hz. 

The  results  of  an  interactive  run  of  program  ACFA  are  shown  in  Fig  2. 

10  PROGRAM  ACFD 


This  program  calculates  the  delay  performance  of  a  given  filter.  The  user 
specifies  the  filter  order  and  filter  type  and  the  program  calculates  the  delay 
at  specified  frequencies  of  interest. 


The  following  input  parameters  are  required  in  the  following  order. 
Unless  stated  otherwise  there  should  be  one  item  per  record. 


(1)  One  of  the  three  symbols  L,  H,  B,  depending  on  whether  the  filter 
is  low-pass,  high-pass  or  band-pass.  The  symbol  should  be  the  first 
character  of  the  record. 


(2)  Filter  name.  This  should  begin  at  the  first  character  of  the  record. 
Only  the  first  four  characters  of  the  name  are  significant,  so  that  a 
Gaussian  filter,  for  example,  could  be  specified  by  either  GAUSSIAN  or 
GAUS. 


(3)  Alpha.  This  real  number  applies  only  to  Ultraspherical  filters 
(see  section  3.5)  and  selects  one  filter  type  out  of  the  Ultraspherical 
family  of  filter  types.  For  other  filter  types  this  parameter  should  be 
omitted. 


(4)  The  order  of  the  filter.  In  the  case  of  band-pass  filters  the  order 
to  be  given  is  the  order  of  the  low-pass  prototype  (ie  half  the  actual 
order) . 
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(5)  The  cut-off  frequency  in  hertz.  This  parameter  should  only  be  given 
if  the  filter  is  low-pass  or  high-pass. 

(6)  Geometric  centre  frequency  and  bandwidth,  both  in  hertz.  These 
numbers  should  appear  on  the  same  record.  The  parameters  should  only  be 
given  if  the  filter  is  band-pass. 

I 

(7)  The  attenuation  in  dB  at  the  cut-off  frequency  (or  frequencies). 

(8)  The  number  of  frequencies  of  interest.  This  is  an  integer  not 
greater  than  20. 

(9)  The  frequencies  of  interest,  in  hertz.  The  numbers  may  be  spread 
over  any  number  of  records. 


10.2  Program  output 


The  delay  at  each  frequency  of  interest  is  output  in  the  form 


DELAY  AT  FREQ.  no.  WITH  ORDER  no.  IS  no. 

The  delay  is  given  in  seconds.  The  order  referred  to  is  the  order  of  the  low- 

pass  prototype  which  in  the  case  of  band-pass  filters  is  half  the  actual  order. 

To  emphasize  this,  in  the  band-pass  case  the  actual  order  of  the  band-pass  filter 
is  output  after  the  last  frequency  of  interest. 

10.3  Restrictions 

Papoulis,  Halpern  and  Synchronously  Tuned  filter  types  are  not  available. 
The  orders  of  the  filter  types  are  restricted  as  shown  in  the  following  table. 

Filter  type  Maximum  order 

Butterworth  30 

Chebyshev  20 

Ultraspherical  20 

Gaussian  8 

Bessel  6 


As  stated  previously  there  can  be  not  more  than  20  frequencies  of  interest. 
If  it  is  desired  to  change  this  limit  to,  say,  50,  then  the  ninth  line  of  program 
ACFD  should  be  changed  from 


[l  :  20]  'REAL'  WINT  ; 
to 

[l  :  50]  'REAL'  WINT  ; 


10.4  Example  of  use 


Consider  a  low-pass  fifth  order  Chebyshev  filter  with  cut-off  frequency 
1  kHz  and  pass-band  attenuation  ripple  of  0.28  dB.  Fig  3  shows  the  results  of  an 
interactive  run  of  program  ACFD,  the  frequencies  of  interest  being  100  Hz,  500  Hz, 
900  Hz,  1  kHz,  1.4  kHz. 

1  I  PROGRAM  ACFC 

This  program  calculates  the  capacitor  and  resistor  values  needed  to  realise 
a  given  filter  type  with  a  given  circuit  configuration.  The  component  values 
are  chosen  where  possible  to  have  certain  desirable  properties,  as  described  in 
section  5.  For  some  circuits  analysis  of  sensitivity  or  circuit  performance  can 
be  carried  out,  as  described  in  sections  6  and  7. 


11.1  Input  parameters 

There  are  many  options  available  to  the  user  of  program  ACFC.  There  is  a 
choice  of  filter  types,  a  choice  of  circuit  types,  a  choice  of  sensitivity  calcu¬ 
lations  and  so  on:  the  data  required  will  depend  on  which  of  the  options  is 
taken.  For  this  reason  the  input  requirements  have  been  illustrated  in  the  flow¬ 
chart  of  Fig  22.  The  flowchart  indicates  input  only  -  it  does  not  show  calcula¬ 
tion  or  output.  Rectangles  indicate  information  to  be  input,  each  rectangle 
corresponding  to  one  record  unless  stated  otherwise  below.  Diamond  shaped  boxes 
refer  to  decisions  which  control  the  path  through  the  flowchart.  Since  some  of 
the  input  consists  of  answers  to  questions  it  is  necessary  to  be  clear  about  the 
distinction  between  a  question  to  be  answered  in  the  data  (rectangle)  and  a 
question  to  be  answered  by  the  person  following  the  flowchart  (diamond) .  A  rect¬ 
angle  containing  a  question  mark  indicates  that  the  data  must  start  with  either 
Y  (for  YES)  or  N  (for  NO) .  A  diamond  shaped  box  containing  only  a  question 
mark  stands  for  the  question  appearing  in  the  immediately  preceding  box.  The 
numbers  appearing  beside  the  boxes  refer  to  the  notes  below. 


(1)  These  are  four  integers.  The  filter  order  is  that  of  the  low-pass 


prototype,  ie  the  actual  order  if  low-pass  or  high-pass  but  half  the  actual 


order  if  band-pass. 


Filter  type 


1  for  a  Butterworth  filter 

2  for  a  Chebyshev  filter 

3  for  a  Gaussian  filter 

4  for  a  Bessel  filter 

5  for  a  Ultraspherical  filter 

0  for  a  Special  filter  defined  by  user 
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Circuit  type 


1  for  Unity  Cain  section 

2  for  Single  Feedback,  section 

3  for  Multiple  Feedback  section 

4  for  State  Variable  section 

5  for  Ring  of  Three  section 

6  for  Mitra-Aatre  section 


LHB 


1  for  low-pass 

2  for  high-pass 

3  for  band-pass. 


(2) 


This  facility 


is  for  use  when  considering  a  filter  which  has  been 


normalised  to  have  unity  delay  at  zero  frequency.  The  user  can  then  de- 


normalise  by  specifying  the  required  delay  at  zero  frequency.  If  the 
filter  is  not  of  this  type  then  the  zero  frequency  delay  should  be  set  to 
unity. 


(3)  The  filter  is  defined  by  the  poles  of  its  transfer  function.  For 
complex  poles  the  real  part  should  be  given,  followed  by  the  imaginary 
part.  The  complex  poles  will  of  course  occur  in  complex  conjugate  pairs: 
only  one  member  of  each  pair  should  be  given. 

(4)  The  order  is  the  order  of  the  low-pass  prototype. 

(5)  Alpha  is  the  parameter  which  specified  one  particular  filter  type 
out  of  the  Ultraspherical  family  of  filter  types. 

(6)  This  is  the  overall  gain  desired  for  the  whole  circuit  and  is  given 
as  a  numerical  value  (not  in  dB) . 


(7)  These  frequencies  should  be  given  in  hertz. 

(8)  The  pass-band  ripple  should  be  given  in  dB. 

(9)  A  convenient  value  of  resistance  is  given  in  ohms. 

(10)  These  are  the  minimum  and  maximum  values  of  capacitance  (in  farads). 

(11)  These  are  the  minimum  and  maximum  values  of  resistance  (in  ohms). 

(12)  The  gain-bandwidth  product  of  the  operational  amplifier  is  given  in 
hertz.  It  is  assumed  to  be  the  same  for  all  the  amplifiers  in  the  circuit. 

(13)  The  overall  gain  has  been  given  earlier,  but  here  the  relative  gain 
of  each  second  order  section  is  given.  (The  first  order  section,  if 
present,  is  not  involved  as  it  always  has  unity  gain.)  The  user  specifies 
the  relative  gains 

V  G2’  G3’  Gm 
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where  m  is  the  number  of  second  order  sections  and  is  equal  to  the  largest 
integer  less  than  or  equal  to  half  the  actial  filter  order.  Gj,  . . . ,  G^ 
must  satisfy 

m 

TK  -  ■  • 

i=  I 

If  the  user  does  not  wish  to  specify  these  relative  gains  they  will  each 
be  set  to  unity. 

(14)  After  the  component  values  have  been  calculated  the  user  may  change 
them  {eg  to  investigate  the  effect  of  using  preferred  value  resistors). 

Once  they  have  been  changed  they  will  remain  changed  and  will  keep  their 
new  values  until  they  are  changed  again.  All  future  calculations,  inclu¬ 
ding  sensitivity  calculations,  will  be  done  with  the  new  values. 

(15)  Some  or  all  of  the  elements  may  be  changed.  Specify  the  number  of 
changes . 

(16)  For  each  element  to  be  changed  one  record  should  be  input  containing 
three  quantities: 

(i)  the  code  identifying  the  element  within  its  section.  This  is 
either  C  (capacitor)  or  R  (resistor)  followed  immediately  by 
the  number  of  the  element  in  the  section, 

(ii)  the  section  number, 

(iii)  the  new  value  for  the  element,  in  ohms  or  farads,  eg 

C2  3  2.3S-9 

R4  1  2000.0  . 

(17)  Answer  YES  if  either  differential  or  large  change  sensitivities  are 
required. 

(18)  Answer  YES  to  obtain  differential  sensitivities  of  each  element, 
considered  separately,  at  two  critical  frequencies.  The  two  critical  fre¬ 
quencies  are  zero  and  cut-off  frequency  (for  low-pass  filters),  infinite 
and  cut-off  frequencies  (for  high-pass  filters),  centre  and  lower  cut-off 
frequencies  (for  band-pass  filters). 

(19)  Give  maximum  percentage  change  in  gain  desired,  so  that  the  program 
may  calculate  the  tolerances  required  to  achieve  this  aim.  If  the  user  is 
uncertain  about  the  value  or  uninterested  in  the  result,  1.0  should  be 
input . 
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(20)  The  simplest  way  to  investigate  sensitivity  is  to  change  some  or  all 
of  the  element  values  and  see  the  change  in  gain  at  given  frequencies.  The 
user  should  specify  the  number  of  elements  to  be  changed.  If  no  elements 
are  to  be  changed,  zero  should  be  input.  It  should  be  noted  that  contrary 
to  the  situation  in  item  14,  the  changes  made  here  are  temporary.  After 
the  large  change  sensitivities  have  been  calculated,  all  elements  revert 

to  their  previous  values. 

(21)  For  each  element  whose  value  is  to  be  changed  a  record  should  be 
input  containing  three  items: 

(i)  the  code  identifying  the  element  within  its  section.  This 
consists  of  either  C  (capacitor)  or  R  (resistor)  followed 
immediately  by  the  number  of  the  element  in  the  section, 

(ii)  the  section  number, 

(iii)  the  percentage  change  to  be  made  in  the  element  value.  A 
percentage  increase  is  assumed  if  the  number  given  is  positive,  a 
percentage  decrease  is  assumed  if  the  number  given  is  negative,  eg 

C2  3  5.0 

R5  1  -1.0  . 

(22)  Answer  YES  to  obtain  the  following  information: 

(i)  input  and  output  impedances  of  each  section,  and  of  the 
complete  circuit, 

(ii)  the  minimum  and  maximum  voltages  and  voltage  spread, 

(iii)  the  transmission  matrix  of  the  complete  circuit, 

(iv)  the  actual  voltage  gain  of  the  circuit  when  using  non-ideal 
operational  amplifiers. 

(23)  For  the  calculation  of  element  values  and  sensitivities  the  ordering 
of  the  sections  is  irrelevant.  For  the  calculation  of  the  information 
referred  to  in  item  (22)  however,  this  ordering  must  be  known.  Answer  YES 
to  this  question  if  the  sections  are  to  be  arranged  in  increasing  order  of 
quality  factor,  Q  . 

(24)  The  order  of  the  sections  should  be  specified,  including  the  first 
order  section,  if  present,  eg 

5  I  3  2  4  . 
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(25)  The  input;  and  output  impedances  of  the  operational  amplifier  should 

be  given  in  ohms  and  the  open-loop  dc  gain  as  a  numerical  value  (not  in  dB) . 
The  output  impedance  should  not  be  given  the  value  zero.  The  same  para¬ 
meters  are  assumed  to  apply  to  all  the  amplifiers. 

(26)  The  information  of  item  (22)  will  be  calculated  at  the  two  critical 
frequencies  defined  in  item  (18).  If  this  information  is  required  at  other 
frequencies  the  number  of  such  extra  frequencies  should  be  input.  If  no 
extra  frequencies  are  desired,  input  zero. 

(27)  Give  extra  frequencies  in  hertz. 

(28)  Answer  YES  to  obtain  either  sensitivity  or  performance  analysis  of 
item  (22). 

(29)  Answer  YES  to  change  either  the  relative  gains  (item  (13))  or  the 
component  values  (item  (14)). 

(30)  Answer  YES  to  consider  a  new  set  of  filter  data. 

11.2  Program  output 

The  program  output  may  be  conveniently  divided  into  four  sections:  transfer 
function,  component  values,  sensitivity,  performance.  The  transfer  function 
information  consists  of  the  poles  of  each  first  and  second  order  second.  For 
each  section  there  is  printed  A  and  B  where  the  denominator  of  the  transfer 
function  of  that  section  is 

2 

s  +  As  +  B  . 

Then  the  quality  factor  Q  ,  and  the  resonant  frequency  wn  ,  are  printed 
(Wq  =  B,  Wq/Q  =  A).  For  a  first  order  section,  a  is  printed  where  the  denomin¬ 
ator  of  the  tt.-.isfer  function  of  the  section  is  s  +  a  .  The  transfer  function 
information  is  given  both  for  the  normalised  transfer  function  of  the  prototype 
low-pass  filter  (under  the  heading  NORMALISED  ROOTS)  and  for  the  de-normal ised 
transfer  function  of  the  actual  filter  (under  the  heading  DE -NORMALISED  ROOTS). 

The  component  value  information  begins  with  the  name  of  the  section  type. 

For  each  section  in  turn  there  is  output  the  results  of  the  iterative  process 
to  find  satisfactory  component  values  for  that  section.  At  the  end  of  each 
iteration  the  program  prints  out  the  values  of  the  capacitances  and  resistances 
generated  in  that  iteration.  If  no  satisfactory  set  of  values  can  be  found  for 
a  section  a  message  is  printed  to  that  effect.  After  the  iteration  results  the 
final  values  found  for  each  section  are  output  under  the  heading  COMPONENT  VALUES, 
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the  values  being  in  farads  and  ohms.  An  indication  is  then  given  of  whether  it 
was  found  possible  to  achieve  the  desired  overall  gain.  The  message  is 
OKGAIN  =  T  if  it  was  possible,  or  OKGAIN  =  F  if  it  was  not  possible.  It 
should  be  noted  that  the  message  OKGAIN  =  F  will  always  be  printed  if  the  unity 
gain  section  is  in  use,  even  if  the  desired  gain  was  in  fact  unity.  If  the 
desired  gain  was  not  achieved,  the  actual  gain  of  each  second  order  section  and 
of  the  complete  circuit  will  be  printed. 

If  differential  sensitivities  were  requested  they  are  printed  out  for  each 
section  at  each  of  the  two  critical  frequencies  (see  item  (18)  of  section  11.1). 
Also  printed  are  the  tolerances  required  for  the  specified  percentage  variation 
in  gain  at  the  critical  frequencies.  Some  elements  have  zero  sensitivities  and 
consequently  infinite  tolerances.  In  this  case  the  tolerance  is  truncated  to 
10*0  x  specified  variation  in  gain.  If  any  large  change  sensitivities  were 
requested  the  resultant  percentage  change  in  gain  at  the  two  critical  frequencies 
are  printed.  A  positive  value  signifies  an  increase  and  a  negative  value 
signifies  a  decrease. 

If  a  performance  analysis  was  requested  (see  item  (22)  of  section  11.1) 
then  the  following  information  is  output.  First,  if  the  sections  were  required 
to  be  ordered  with  increasing  quality  factor  then  this  order  is  printed.  Then 
for  the  two  critical  frequencies  and  for  any  extra  frequencies  requested,  the 
program  outputs  the  absolute  values  of  the  input  and  output  impedances  of  each 
section,  the  minimum  and  maximum  voltage  and  voltage  spread,  the  transmission 
matrix  of  the  complete  circuit,  the  modulus  and  phase  of  the  input  and  output 
impedances  of  the  complete  circuit  and  the  modulus  and  phase  of  the  voltage  gain 
of  the  complete  circuit. 

Whenever  a  section  number  is  referred  to  in  the  output  (or  the  input)  it 
generally  refers  to  the  number  as  printed  in  the  component  values  information 
under  the  heading  COMPONENT  VALUES.  There  is  one  exception  to  this  rule.  The 
maximum  and  minimum  voltages  are  given  in  the  performance  information  by  messages 
of  the  form 

MINIMUM  VOLTAGE  (AT  END  OF  SECTION  nn)  IS  xxx 

MAXIMUM  VOLTAGE  (AT  END  OF  SECTION  nn)  IS  xxx  . 

Here  nn  refers  to  the  number  relative  to  the  current  ordering.  For  example,  if 
there  are  four  sections  which  have  been  put  in  the  order 
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then  "END  OF  SECTION  2"  refers  to  the  end  of  the  second  section  in  the  new  order¬ 
ing,  ie  that  originally  numbered  1.  Note  that  "END  OF  SECTION  0"  is  printed 
instead  of  "START  OF  SECTION  1". 


As  mentioned  in  section  8,  the  output  will  also  contain  the  questions  asked 
by  the  program  to  elicit  the  data  from  the  user  in  interactive  mode.  To  remove 
these  questions  from  an  output  file  the  macro  ACFTIDY  can  be  used. 


11.3  Restrictions 


Only  the  following  filter  types  are  available:  Butterworth,  Chebyshev, 
Ultraspherical ,  Gaussian,  Bessel,  User-defined.  The  orders  available  for  each 
filter  type  are  restricted  as  indicated  in  the  following  table. 


Filter  type 

Butterworth 

Chebyshev 

Ultraspherical 

Gaussian 

Bessel 

User-defined 


Maximum  order 
30 
20 
20 
8 
6 

20  . 


Any  of  the  circuit  types  of  section  4  can  be  chosen  though  it  should  be  noted 
that  not  all  of  the  types  are  available  in  low-pass,  high -pass  and  band -pass  form. 
The  following  table  shows  which  combinations  are  available 


Low-pass  High-pass 
Unity  Gain  /  / 

Single  Feedback  /  / 

Multiple  Feedback  /  / 

State  Variable  /  / 

Ring  of  Three  / 

Mitra-Aatre 


Band -pass 

/ 

/ 

/ 

/ 

/ 


Sensitivity  and  performance  information  is  available  only  for  Unity  Gain 
and  Ring  of  Three  circuit  types. 

The  standard  capacitor  values  used  are  listed  in  section  5.  It  may  be 
desired  to  change  these  values.  Program  ACFC  can  be  edited  to  accomplish  this, 
in  the  following  way.  If  more  than  20  standard  values  are  required,  then  the 
eleventh  line  of  the  program  should  be  changed  from 
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to 


[l  :  20]  'REAL'  CASTAN  ; 
[1  :  nn]  'REAL'  CASTAN  ; 


where  nn  is  the  number  of  standard  values  required.  The  following  seven 
consecutive  lines  control  the  setting  up  of  the  standard  values  in  the  array 
CASTAN. 

NCASTAN  :=13  ; 

'FOR’  I  ’TO’  4  ’DO’  ’BEGIN’  'FOR'  J  'TO'  3  'DO'  'BEGIN' 

'REAL'  TM1 ,  TM2  ;  TM1  :=  'CASE'  I  'IN'  1.0&-11,  1.0&-10,  1.0&-9,  1 .0&-8 

'ESAC'  ; 

TM2  :=  'CASE'  J  'IN'  1.0,  2.2,  5.0  'ESAC'  ; 

CASTAN [3* ( I- 1)  +  J]  :=  TM1  *  TM2 

'END'  'END'  ; 

CASTAN  [l  3]  :=  1.0&-7  . 

These  lines  may  be  changed  or  replaced  as  required. 

11.4  Examples  of  use 

Consider  a  twelfth  order  band-pass  Bessel  filter  with  centre  frequency 
500  Hz  and  bandwidth  20  Hz.  The  overall  gain  is  desired  to  be  unity  at  the 
centre  frequency.  A  Ring  of  Three  implementation  is  required  with  capacitances 
between  10  pF  and  0. 1  pF  and  resistances  between  1  kft  and  1  MO.  A  convenient 
arbitrary  resistance  value  is  10  kfl.  Differential  sensitivities  are  required, 
with  tolerances  to  be  calculated  for  a  maximum  10%  change  in  gain.  Fig  4 
illustrates  a  non-interactive  run  of  program  ACFC,  the  data  being  shown  in  Fig  5, 
and  the  results  in  Fig  6.  It  can  be  seen  that  four  iterations  per  section  are 
required  to  find  satisfactory  component  values,  the  arbitrary  capacitance  being 
reduced  step  by  step  until  all  the  resistances  are  less  than  1  MO.  Two  of  the 
resistances  in  each  section  have  the  given  arbitrary  value  of  10  kO,  but  the 
other  resistances  will  have  to  be  rounded  to  convenient  values.  It  may  be  anti¬ 
cipated  from  the  differential  sensitivities  that  this  rounding  will  have  most 
effect  in  sections  3  and  6,  the  resistors  R^ ,  R2,  R^,  R^  being  particularly 
sensitive. 


Fig  7  shows  an  interactive  run  of  program  ACFC.  In  this  case  we  are  con¬ 
cerned  with  a  fourth  order  Butterworth  low-pass  filter  of  cut-off  frequency 
200  Hz.  Again  a  Ring  of  Three  implementation  is  required  with  resistance  and 
capacitance  bounds  as  in  the  previous  example.  The  component  values  are  found 
in  one  iteration  for  the  first  section  and  two  iterations  for  the  second  section. 


42 


On  this  occasion  we  do  not  ask  for  differential  sensitivities,  but  we  investigate 
the  effect  of  using  5%  resistors  and  2%  capacitors  for  the  first  section.  We  find 
that  we  might  obtain  a  decrease  in  gain  of  almost  5%.  Assuming  typical  values  of 
I  Mil  for  the  input  impedance,  300  0  for  the  output  impedance  and  100  dB  for  the 
open-loop  gain  of  the  operational  amplifier  we  find  that  the  input  impedance  of 
the  complete  circuit  is  about  800  kfi  and  the  output  impedance  is  of  the  order  of 
milliohms.  The  gains  at  zero  and  cut-off  frequency  are  approximately  0  dB  and 
-3  dB,  as  required.  We  then  decide  to  do  another  large  change  sensitivity 
calculation,  this  time  involving  only  the  second  section.  We  see  that  a  decrease 
in  gain  of  2%  is  predicted.  Of  course,  before  concluding  that  the  second  section 
is  less  sensitive  than  the  first  we  should  try  changing  the  components  in  differ¬ 
ent  senses  and  also  look  at  the  differential  sensitivities.  But  we  decide  that 
we  have  enough  information,  and  terminate  the  run. 

12  CONCLUSIONS 

We  have  described  haw  a  low-pass,  high-pass  or  band-pass  filter  may  be 
designed  by  first  using  program  ACFA  to  ensure  satisfactory  attenuation  perform¬ 
ance,  then  checking  the  delay  performance,  if  necessary,  using  program  ACFD,  and 
lastly  obtaining  component  values  from  program  ACFC.  For  some  circuit  types, 
program  ACFC  can  then  be  used  to  investigate  the  sensitivity  of  individual 
components  or  groups  of  components  and  check  the  circuit  performance  with  non¬ 
ideal  operational  amplifiers.  We  have  described  how  these  further  calculations 
are  performed  for  the  circuit  types  in  question  and  also  how  they  may  be  performed 
for  other  circuit  types. 
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Appendix  A 

GAIN  AND  FREQUENCY  TRANSFORMATIONS 
(see  section  2) 

To  change  the  gain  from  unity  to  G  ,  say,  replace  the  term 


Ck8  +  Dks  +  % 

S'  +  V  +  \ 


W  +  W  +  GkEk 
s2  +  V  +  \ 


where  Gj,  G^t  G^  ...  are  such  that 


11/  L. 

tt*- 


if  n  is  even,  or 


(n-l)/2 


°k  " 


if  r<  is  odd. 


A. 2  To  generate  a  low-pass  filter  of  cut-off  frequency  W  radians  per  second 
replace  s  by  s/W  ,  ie  replace 


_e _  ,  eW 

+  a  ^  s  +  aW 


and  replace 


3  +  A^  +  ^ 


S2  +  AkWs  +  ^w2 


A. 3  To  generate  a  high-pass  filter  of  cut-off  frequency  W  radians  per  second 
replace  s  by  W/s,  ie  replace 


s  +  a 


es/a 
s  +  W/a 


and  replace 
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+  V  + 


by 


E,  /  B, 
k  k 


s2  +  A^Ws/Bk  +  V^/B 


A. 4  To  generate  a  band-pass  filter  of  lower  and  upper  cut-off  frequencies 
and  Wg  replace  s  by 

»2  ♦  V, 

<"b  -  Vs  ' 

This  will  generate  a  transfer  function  which  has  order  double  that  of  the  low- 
pass  prototype  and  is  geometrically  symmetric  about  its  centre  frequency 
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ALGORITHM  FOR  SATISFACTORY  COMPONENT  VALUES 
(see  section  5) 

From  section  4.6  the  design  equations  for  the  Ring  of  Three  low-pass 
circuit  are 


0)  =  R2/(R]R3R4c)c2) 


Q  =  R5C2co0 


G  =  R4/R6  . 


Having  set  R^ ,  R2  and  Cj  =  C2  we  calculate 


Cl“WR| 


Q/(c,o>0) 


r3/g 


R4  =  R3 


00 

o 


Suppose  we  find  that  R«  >  R  .  We  ask:  can  we  reduce  R~  by  changing  one  of 

•j  max  j 

the  arbitrary  resistances?  The  answer  is  that  we  can,  by  suitably  re-setting 

R,  .  In  order  not  to  change  R„  too  much  from  its  previous  value  we  set 

Z  Z  .222 
R0  =  R  and  calculate  the  corresponding  value  of  R~  ,  which  is  C.to-.R  R,  . 
3  max  re  2  10  max  1 

Now  suppose  instead  that  R^  >  Rmax  •  Again  we  can  change  R^  by  changing 

R„  ,  and  the  value  of  R„  which  makes  R,  =  R  is  C^u^R^  GR,  .  To  allow 
2  26  max  1  0  max  1  222 

for  the  possibility  that  R_  >  R  and  R,  >  R  we  set  R.  to  C.w„x  R. 

3  max  6  max  2  10  1 

where  x  =  min(R  ,  R  G)  .  Similarly  if  R,  <  R  ,  or  R,  <  R  .  . 

max  max  3  mm  6  min  we  calcu¬ 


late  the  value  of  R„  necessary  to  make  R„  =  R  .  or  R, 

2  3  mm  6 


=  R  .  . 
mm 


Now  suppose  that  R^  > 


max 


In  this  case  we  cannot  change  R^  without 


changing  Cj  .  From  the  design  equations  we  see  that  in  order  to  reduce  R^  we 
must  increase  Cj  .  This  is  not  the  only  reason  for  which  we  might  want  to 
increase  C 


'1  * 


If  R2  has  been  reduced  because  R3  or  R^  were  out  of  bounds 


as  described  above,  then  it  might  be  that  the  new  value  of  R2  was  less  than 

R  .  .  Then  C.  should  be  increased  so  that  in  the  next  iteration  R~  and  R, 
min  1  3  6 

will  be  decreased.  Also,  we  must  check  that  C.  has  not  become  less  than  C 

1  max . 
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Thus  we  must  increase  to  the  next  higher  preferred  value  if 

R,  >  R  or  R„  <  R  .  or  C.  <  C  , 

5  max  2  mm  1  mm 

Similarly  we  must  decrease  if 

Rc  <  R  .  or  R„  >  R  or  C,  >  C 

5  mm  2  max  1  max 

If  Cj  need  not  be  changed  then  we  have  succeeded  in  calculating  a  satisfactory 
set  of  component  values  for  this  section.  If  Cj  is  changed,  then  we  re-set 
Rj  and  R2  to  the  preferred  value  and  repeat  the  procedure.  We  can  fail  to 
obtain  a  satisfactory  set  of  component  values  for  any  of  the  following  reasons. 

(i)  We  may  be  unable  to  increase  C.  because  that  would  make  C,  >  C 

1  1  max 

(or  unable  to  reduce  it  because  that  would  make  C,  <  C  .  )  . 

1  min 

(ii)  Our  algorithm  may  tell  us  to  increase  and  decrease  Cj  simultaneously 

eg  if  we  have  both  Rc  >  R  and  R„  >  R 

5  max  2  max 

(iii)  Suppose  at  a  certain  iteration  we  decide  to  increase  Cj  and  at  the 
next  iteration  we  decide  to  decrease  Cj  .  Then,  unless  the  situa¬ 
tion  were  detected,  the  algorithm  would  oscillate  forever  between 
two  preferred  capacitance  values. 
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For  a  low-pass  filter  F^(s)  can  written  38 


Gu! 


so 


Whence 


Fl(s) 


3Fl(s) 


dp 


2  /  2 
s  +  u>0  s/Q  +  uQ 

G|(s2  +  ioqs/Q  +  uQ^2a)Q  ^ 

-4q£-^V4 

+  2“o  3r|/(s2  +  “oS/Q  *  0 

+  “o  lf/(s2  +  v/Q  +  “o) 


)j(s2  +  UqsJ Q  +  co2)2  ^./<o 


K 


q^£-w  H 

w  3p  0  3p 


)/q2]s  +  2/wo  3irj fa2  +  W0s/Q  +  “o, 


Similar  formulae  can  be  obtained  for  other  types  of  filter. 


( 

I 


G  9p 


038 


038 
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DERIVATION  OF  TRANSMISSION  MATRIX 
(see  section  7) 

With  non-ideal  operational  amplifiers  the  second  order  Ring-of-Three  section 
is  modelled  as  shown  in  Fig  21.  Considering  the  nodes  from  2  to  10  in  turn  we 
have  the  equations 


where  R^n>  Rout’  are  t*ie  *nPut  impedance,  output  impedance  and  open-loop 
gain  of  the  amplifier!  and  comparing  Figs  19  and  21  we  see  that 


Z,  =  R5/0  +  R5C2s) 

Z2  =  1/CCjS)  . 
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Manipulation  of  these  equations  gives 


Call  these  equations  (7-1).  We  wish  to  determine  the  elements  of  the  open 
circuit  impedance  matrix 


Z 

'Z.l 

Z 1 2 

1  ’ 

04 

ESI 

Z22 

1 

where  Z  j  = 

V, 

9 

Z12 

= 

h 

9 

I. 

m 

Iout=0 

*out 

Iin=0 

II 

04 

N 

V10 

9 

Z22 

= 

vjo 

• 

I. 

in 

^"out=0 

■'"out 

*in=0 

(The  positive 

inwards  current  convention  is 

used . ) 

We  have 

I. 

in 

(V! 

1  -  V'R6  • 
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5) 


De  f ine 

“i  ■  Vv, 

“2  '  VV, 

“3  '  VV. 

“4  '  VVI 

“5  '  VV1 

“6  ‘  V.0/Vl  ' 

Consider  first  the  situation  when  I  =  0  .  We  can  re-write  equations  (7~1) 

out 

in  the  form 

Mu  =  b 


where  M  is  a  6  *  6  matrix  and 


]/r6  ,  b2  =  b3  =  b4 


b5  "  b6 


These  equations  can  be  solved  for  u  and  we  shall  then  have 


zn  ■  viVlvi  -  V  ■  V'1 


z2l  ■  <vio'V/(W  ■  "6zn 


Now  consider  the  situation  when  I^n  =  0  .  Equations  (7-1)  will  still  apply, 
except  that  we  should  replace  the  last  equation  by 


V  =  V 
1  2 


z-e  u„ 


So  we  Still  have  Mu  =  b  with  M  being  unchanged  except  for  its  last  row  and  J) 
being  unchanged  except  that  now  b,  =  1  .  From  the  last  equation  of  (7-1) 

D 


i  /v. 

out  1 


whence 


Z,,  =  V,/1  _ 

12  1  out 


and 


J22 


U6Z12 


can  be  found  . 
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The  transmission  matrix  can  now  be  found  from  the  relations 


all  “ 

zn/z2> 

ai2 

'  1/Z2I 

a2 1  " 

where  *  -  Z,  ,2^  -  • 

CM 

SJ 

a22 

'  Z22/Z 

o 

UJ 
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LIST  OF  SYMBOLS  AND  ABBREVIATIONS 


al 1 ,ai2’a21 ,a22 
A 

A0 

\ 

ACFA 

ACFC 

ACFD 

B 

Bk 

Bn(x) 

C 

Ck 

C j ,0^  etc 
C 

max 

C  . 
min 

d 

D 

°k 

D(w) 

e 

E 


negative  of  pole  of  transfer  function  of  first  order  section 
elements  of  matrix 

coefficient  of  s  in  denominator  of  second  order  transfer 
function 

operational  amplifier  open-loop  gain 

coefficient  of  s  in  denominator  of  transfer  function  of 
kth  second  order  section 

transmission  matrix 

name  of  program  to  calculate  attenuation  performance  of 
active  filter 

name  of  program  to  calculate  component  values  of  active  filter 

name  of  program  to  calculate  delay  performance  of  active 
filter 

constant  coefficient  in  denominator  second  order  transfer 
function 

constant  coefficient  in  denominator  of  transfer  function 
of  kth  second  order  section 

nth  order  Bessel  polynomial 
2  . 

coefficient  of  s  in  numerator  of  second  order  transfer 
function 

2  . 

coefficient  of  s  in  numerator  of  transfer  function  of  kth 
second  order  section 

capacitance  values 

maximum  capacitance  value 

minimum  capacitance  value 

coefficient  of  s  in  numerator  of  first  order  transfer 
function 

coefficient  of  s  in  numerator  of  second  order  transfer 
function 

coefficient  of  s  in  numerator  of  transfer  function  of  kth 
second  order  section 

delay  at  frequency  w 

constant  coefficient  in  numerator  of  first  order  transfer 
function 

constant  coefficient  in  numerator  of  second  order  transfer 
function 

constant  coefficient  in  numerator  of  transfer  function  of 
kth  second  order  section 


V> 

G 

H(s) 


transfer  function  of  kth  second  order  section 
absolute  value  of  gain  of  filter  or  filter  section 
filter  transfer  function 


j 


54 


Rj,!^  etc 


T„« 

Uj,u2,  etc 
"»<»> 

Vj,V2,  etc 
Y^(x) 

ZH,Z12’Z21,Z22 


zrz2 


LIST  OF  SYMBOLS  AND  ABBREVIATIONS  (concluded) 

input  current 
output  current 
filter  order 
arbitrary  parameter 

operational  amplifier  gain-bandwidth  product 
quality  factor 

variable  measuring  the  performance  of  a  circuit 

resistance  values 

maximum  resistance  value 

minimum  resistance  value 

Laplace  transform  variable  =  jw 

differential  sensitivity  of  variable  r  with  respect  to 
change  in  parameter  p 

nth  order  Papoulis  polynomial 

voltage  ratios 

nth  order  Halpern  polynomial 
voltages 

nth  order  Ultraspherical  polynomial 

elements  of  Z  matrix 

open  circuit  impedance  matrix 

impedances  used  in  calculating  non-ideal  performance  of 
Ring  of  Three  circuit 

parameter  relating  to  Ultraspherical  filters 
attenuation  at  w  =  1 

time  constant  used  in  the  definition  of  State  Variable 
design  equations 

arbitrary  frequency  in  radians  per  second 

2 

resonant  frequency:  =  B 
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«-  COMPILE  ACFA 


Fig  t 


«-  EXECUTE 
•  «  • 

TYPE  L,  H  OR  B  FOR  LOW,  HIGH  OR  BAND-PASS 
«-  H 

TYPE  FILTER  NAME 
BESS 

TYPE  ORDER  OF  FILTER 
4 

TYPE  FREQUENCY  OF  PASS  BAND  EDGE 
«-  1.5&3 

HOW  MANY  FREQUENCIES  OF  INTEREST? 
i-  3 

TYPE  FREQUENCIES  OF  INTEREST 
«-  1.2&3  1.0&3  0.8&3 


ATTENUATION  AT  FREQ. 

1 .20000000&3 

WITH 

ORDER 

4 

IS 

4 . 9398 1 97 1 39&0 

ATTENUATION  AT  FREQ. 

1 . 00000000&3 

WITH 

ORDER 

4 

IS 

7 .421822887040 

ATTENUATION  AT  FREQ. 

8.00000000&2 

WITH 

ORDER 

4 

IS 

1 . 1 84928468 1  &  1 

Fig  1  Example  1  of  ACFA 


COMPILE  AC FA 


•e  EXECUTE 

TYPE  L,  H  OR  B  FOR  LOW,  HIGH  OR  BAND-PASS 
B 

TYPE  FILTER  NAME 
CHEB 

TYPE  ORDER  OF  FILTER 
+-  0 

TYPE  ATTENUATION  AT  STOP  BAND  FREQUENCY 
f-  70.0 

TYPE  CENTRE  FREQUENCY  AND  BANDWIDTH 
+-  500.0  50.0 

TYPE  ATTENUATION  AT  PASS  BAND  EDGE 

+-  0.  1 

HOW  MANY  FREQUENCIES  OF  INTEREST? 

■*-  1 

TYPE  FREQUENCIES  OF  INTEREST 
+■  400.0 

ATTENUATION  AT  FREQ.  4.0000000042  WITH  ORDER  1  IS  1.678143008340 

BAND-PASS  ORDER  IS  +2 

ATTENUATION  AT  FREQ.  4.00000000&2  WITH  ORDER  2  IS  1.572207939341 

BAND-PASS  ORDER  IS  +4 

ATTENUATION  AT  FREQ.  4.0000000042  WITH  ORDER  3  IS  3.457990826141 

BAND-PASS  ORDER  IS  +6 

ATTENUATION  AT  FREQ.  4.0000000042  WITH  ORDER  4  IS  5.355397178641 

BAND-PASS  ORDER  IS  +8 

ATTENUATION  AT  FREQ.  4.0000000042  WITH  ORDER  5  IS  7.252952752141 

BAND-PASS  ORDER  IS  +10 

Fig  2  Example  2  of  ACFA 
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•*-  COMPILE  ACFD 


EXECUTE 

TYPE  L,  H  OR  B  FOR  LOW,  HIGH  OR  BAND-PASS 
•<-  L 

TYPE  FILTER  NAME 
■*-  CHEB 

TYPE  ORDER  OF  FILTER 
5 

TYPE  FREQUENCY  OF  PASS  BAND  EDGE 
1.0&3 

TYPE  ATTENUATION  AT  PASS  BAND  EDGE 
+  0.28 

HOW  MANY  FREQUENCIES  OF  INTEREST? 

5 

TYPE  FREQUENCIES  OF  INTEREST 

100.0  500.0  900.0  1.0&3  1 .48.3 

CHEBYCHEV  FILTER 


DELAY 

AT 

FREQ. 

1 .000000008.2 

WITH 

ORDER 

5 

IS 

6.1) 063896 16&-4 

DELAY 

AT 

FREQ. 

5.000000008.2 

WITH 

ORDER 

5 

IS 

6 . 7223884753&-4 

DELAY 

AT 

FREQ. 

9.00000000&2 

WITH 

ORDER 

5 

IS 

9 . 86367282238-4 

DELAY 

AT 

FREQ. 

1 . 000000008.3 

WITH 

ORDER 

5 

IS 

1 .44 1 48733398-3 

DELAY 

AT 

FREQ. 

1 . 400000008.3 

WITH 

ORDER 

5 

IS 

2 . 63872884938-4 

Fig  3  Example  of  ACFD 


*■  COMPILE  At  FC 


«-  GETEBM  PROGC 
•*-  AS  *TRO,  LJRDAT 
f  AS  *LPO,  LJRRES 
<-  EXECUTE 

ACFTIDY  LJRRES 


Fig  4  Non-interactive  run  of  ACFC 


I^SSm 
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«-  LF  LJRDAT 


6  4  5  3 

1.0 

500.0  20.0 

1.0&4 

1.04-11  1.04-7 

1.043  1.046 

N 
N 

Y 

Y 

10.0 

0 

N 


N 


N 

N 


Fig  6  Data  for  non-intaractiva  run  of  ACFC 
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Fig  8  Results  of  non-interactive  run  of  ACFC 


*k«.  1 1  Oft  ? 

C  AP81'  *  TANCFS  1  .  lit.  n:nii,  -J»  1  .(iiT0’>O‘t*  -rt 
RESlST*Nl|  S  1.MI.U  IIIIIN  44  1  .  nOofluO*  *4  3. 

section  < 

CAPACITANCES  1,00  >000*  I.OOftOO'i*  -V 

SfSISTA'.CK  1  ..il-tjoiiiiF  44  I.OCi.iOOOn  *4  3. 

SECT  3 

CAPACMwNftb  ?,2"  OOP  -C  -V 

RESISTANCES  1  .  OO'J.'Hu*  *4  1  .  (lOnOli  Ur.  *  4  1. 

SECTION  « 

CAPAC  »  t  a  w  i;  f  J*  %  ,  oonooop  -o  S.oomiloi*  -9 
RES  I  ST/»N(.|.S  I.Onh.iftuA  ♦*  I.imW.Ooo*;  *4  to. 
SCrT  !•»»:  3 

CAPAC  I  TANf.i.  J>  -H  1.0f*UW«uA  -8 

RESISTANCES  1.  non  .110*  ♦*»  1,oi.oni*o,  *4  3. 
SFCTIliM  4 

CAPACITANCES  1  ,UU"OU(K  -9  l.fo'Ouuu*  -V 

RFSISTamOS  I.iomi.iOoa  ♦*.  1  .  OfiOOuO*  ♦A  3, 

SECTION  4 

CAPAC  ITANCf!»  -4  ?.?,ino<MIS  -4 

resistances  i  *a  l.onnooui.  *4  1. 

SECTION  4 

-  CAPAf 1 1 AUCFS  S.i»u.M»yu>,  -9  S.O.inouti*  -V 
RESISTANCES  1.000  *00*  l.tiiViflOUN  *4  6. 
SECTION  4 

CAPACITANCES  1  JOi.:T'iMJa  -8  I.CmOoiM*  -8 
RESISTANCES  l.itOO  iuu*  *4  I.UPoOuOr.  *4  3. 

SECTION  s 

CAPAC  ITAHCFS  I.UO.'OiMii.  »9  I.AflAuUo'i  -V 
RESISTANCES  1.0U0-m«A  44  I.OOOOOOu  44  3. 
SECTION  s 

Capacitances  x.«:o>'No<i^  -o  ?.?.inuiin*  -o 
RES! STAMtl S  1.000. HOIS  *4  I.UOoOOok  44  1. 
SFCT10N  S 

CAPmCI  f  A  N  C  F  5  «,  .Oli-MIOM.,  -0  S.PonulMiH  "v 
RESISTANCES  I.OUO'IHIA.  44  I.oO.iIWU*  44  to  . 
SECTION  S 

CAPACITANCES  1  .  OUiinmi,  -8  I.OimVMim*  -8 
RESISTANCES  I.IIUll.lOON  44  1 . 00»i  OOI'N  44  i. 
SECTION  A 

Capacitances  1  ,uo  »<»u"r,  ~v  i.o.hiouo*  -o 
RESISTANCES  1.00IHI0US  44  1.UO«.f»oO<,  44  3. 
SECTION  6 

CAPACITANCES  i’./O/MlOO#:  -4  ?.?OI»IU>H<  -O 
RESISTANCES  1.ii00»90ft  44  I.OOodOOu  *4  1. 
SECTION  A 

CAPAC  T  TAfli.  FS  S.OIMOO"*  -«  S.OOftOIOIA  -O 
RESISTANCES  1.000.01014  44  1.0ftolidU\  44  6. 
SECTION  f> 

CAPAC  IT  ANlFS  1  ,  OUiifMHUi  -ft  l.oOOooi'S  -8 
RESISTANCES  1  .  OOIH'OUA.  44  1  .  OiloOOOA  44  3. 

COMPONENT  V  a  tOE  S 


44  3.?4557nn 


♦5  3 .2906734 


4  S  1.495760* 


44  6  (  S(t1  3^6A 


44  3.290673* 


♦ 5  J. 162725* 


.♦5  1.43740?* 


♦4  6.325*51* 


♦4  3. 1A272S* 


♦  S  3, 1216. SIS 


♦  S  1.41V014* 


♦  4  6.243661* 


44  3 . 1 ?1 8*  1  X 


♦S  3.079041* 


♦s  1.3905444 


♦4  6.158083* 


♦  4  3 . 079041  A 


♦  5  4.803503*  4.5 


46  4,319132*  *6 


♦  a  1  .063242*  46 


46  A , 6 3 8263*  45 


43  4.3101  32*  •  3 


46  4.029683*  *6 


46  2.240763*  46 


♦6  9.85QJ67*  45 


45  4,929683*  *S 


>6  4.620366*  46 


46  2.100166*  46 


46  9.240732*  45 


45  4.620366*  4« 


46  4.041357*  46 


46  1 . B36980*  46 


46  8.082714*  45 


4 S  4.041  357*  45 


SECTION  1 

CAPAC  TTANCfS  1  ,Ul»TliiiH*  -8  1.( 

RESISTANCES  1.O0W.MM1*  *4  1.01 

SECTION  / 

CAPAC  1 T ANTES  1,001000.,  -n  1.1 
RES  IS  T  A  N  C  |  S  1  .  liUUoOON  *4  1.0(i 
SECTION  3 

CAPAC  I1ANIEN  1  .UOmIIOo,,  -8  1.1 

RESIST*. < irs  1,'iuo  i.i(D.  4 it  i.u( 

SECTION  4 

CAPACITANCES  I.UOiOOO.  -6  1.1 
resistances  i,(nt'H>iws  *t.  i.m 


11O000*  -0 

■  1 II  011  4  44  3.2MoUTx 


itOUtio*  -  rt 

1 1 110 ik.  *4  3. 2*. 557c* 


nOoo.i*  -h 

01*011.,  +  4  3  .  0o7  3* 


HOtJOUX  -o 

"iH'W.i  *4  5.1o272S* 


♦•w  3.203603*  a/.  .  S  >  *  *  h  7i  -.OO’tOUS  i 

*u  3.2«,S57 04  44  5.."72^sv*  45  4.803503*  <5 

44  3..>v«*6’3.<  -4  8.644540*  45  4.319132*  ♦! 

♦  «.  3. 1.6?  7  ?5  a  *4  5.031^14*  ♦  b  ^.9:968.11  *5 


Fig  6  contd 


I 


Fig  6  contd 


S ‘I  '  I  1  ■ 1 

Cl,--  ,'■■■■  ■ .  .  . . 

‘  *■  1  '  l  -  ■  •  -  1  . . 

SFf  T 

CM>’  >  »  '  iT»  1  •  -  *  * . .  *»» 

M  i  *,  1,i.  t  .  •>*!■.  Ii  .i  *4  J.M'O****# 

0,r.^r  -  > 


♦  4  1,1/lMlft  *u  •i.a-ams/*  **>  4.h/M6d*  *5 

♦  -  1.17VM41*  *4  ♦*»  4,0*11*/*  *5 


0|**»  »■ • 
**.»►■»  t  i 
S  •>  r  F I »  ' 

C  t 

c  < 

N  1 

a  / 
a  * 
a  <• 
a  > 

a  ft 

S€»  ri'*\ 

c  i 

c  / 

a  i 
a  / 
a  * 
a  * 
a  > 
a  & 

S*r  r 
C  1 
C  / 
a  i 
a  / 
a  3 
a  <. 
a  S 
S  o 
s  e  r  r  r  <  i  r: 
C  1 
C  / 
a  1 
a  * 
a  3 
a  4 
a  5 
a  t 

Sfrrjoft 

c  i 

C  i 

a  1 
a  / 
fl  3 
a  4 
a  5 
a  6 
SICTin* 
C  1 

c  < 

a  1 
a  / 
a  3 
a  4 
a  ) 
a  ft 


I  Ml  -  ’  i'.  ;  I  I  •  T  |l 
.1  ,  .1  I  l  *»l  I  / 


-  \  .  1  •.•<•»»  •  i/«. 

-1.1..  ••■'  •  •  •" 
-1  .1  I-*"  ./  .  *  ■ 

1  .1  .  '•  </  «• 

-  1.  1"  •  .•  •/•  *•' 

-  <  .  1  i  •  %  . 

W  .  V  •••-.  .«)  >  -1 

/ 

-ft.  3r»1/i 

1  .V  -14*.  ♦>' 
-ft  .3'  1/1 

H. ISVIfi/f. 

-ft.  3*  1/1  .»./•.  *•. 
-ft.l'.I.M  >•/’.  *» 

ft  .  ft  /'.V  »■>  ^  -1 

-1  )>■•»«•  .«»>•«.  *o 

-1.1  .«/•/*  •/*•».  ♦  I 

-1  .  1  .14  I  l«  ./ 4.  ♦! 

-1  . 1.,/Vi.  »/i>*  ♦! 

I. Iu’hV"'  *1 
-1  .  lil/Vfc’VU..  *1 

-1.1  ii  /V 4  //«•»,  *1 

?  .5?4/i.  »sv  v  -1 
-1  .  lift. him 
4 

1.  14"-'V  •Ur*.  ♦  <• 

1 . 1  (IHI*/,'4‘>»  ♦" 

I.U'MiVI'ft  *»' 

—  1.1  4<M»V  ’  1  Uft  *»* 

1 .  1  4ft'  v  .*1  u», 

1 . 1  4 *1  u*  *n 
‘i.V-v/liU  -I 
-1  .  M«0*.  •!» 

fl.  A.1. ///«./ 3  A  *V 

H. JM.'I  •/*>.  ♦  .» 

H  .ft.-///-/}* 

«.A.‘V//-/V 
ft  .  A. I.*//-/*  . 
ri . ft,.1"  Aft  1 4ft  a  -  1 
- 1  .  ortili'u  n»*i*.  ♦)> 

1  *1 
1  .  1  I'lfti  .  *>1 »  -1 

1  .  1,. -Ill  .M/ft  ♦'! 
-1.1  4  1  "  .M/A  *  1 

I. 1-1-  >1/  >  -I 

1.1  4  1  •  r-*/  • 

/  .  3.  -/*.  .v.jft  -1 


AT  r,  ftfn  F  *»«»•»*  ►TV 
T'>U.«AM.»  <  X  1 


3  .  /  >.  ♦  u 
3.  /ft 


3 . / ft  *0 
l./ft  *i» 
1.0ft  ♦ 1 

1  . uft  *1 

l./ft  ♦<> 

l./ft  *v 

1  ./ft  ♦!/ 

1.<M  *'l 

1  .^4  *0 

1  ./ft  ♦'! 
1  . ft  ft  *1 
1  . Mft  *1 

V.Hft  -1 

fl.44  -1 

V.U*  -1 

V.ftfc  -1 
V.Mft  -1 
V.ilA  -1 
4  .  /  4 
1.i»S  *1 

i./A  ♦() 
5. /ft  *0 
l.iH  *rj 
3. /A  *0 
3./S  41) 
3. /ft  ♦«! 

1  .lift  *1 
1.*'»  ♦! 

l./ft  ♦.) 

1  ./ft 
l./ft  *i> 
1  ./A  40 
l./ft  *0 
1  ./a  ♦() 
1  .ftft  *1 
1  .  i*ft  *1 

n  .  4  >c  - 1 
V  .  lift  -  1 
ft . 4ft  -1 
ft  .  4ft  -1 
ft  .  4ft  -  1 
r  . 4ft  -i 
4  .  /ft  *1 
1 .  II  ft  *  1 


I 


Fig  6  contd 


Fig  6  concld 


ii»i  "  i.\  i  r  i  v  * 

t  |l< 

A  T  C‘,,,»F* 

)N  *  1 1 1  •  f  ft  c  v 

lf*MI 

-►  .s  r  r  i  v  i  t  • 

Ti.lM.ANLk  (  X  1 

CCTI.Mi 

1 

i 

S.*21  in 

♦  II 

1  .  /ft 

♦U 

t 

S  .  /  I  A- 

•  > 

3  ,  ft* 

♦  1) 

1 

iln- 

... 

i .  2ft 

♦  0 

2 

->.«  *1-  i. 

i .  t>< 

♦  ll 

s 

4  .  C  '1 1  n  .so*. 

♦  .' 

y .  2* 

♦  ll 

4 

S  ,  n"*1  on  ,*»hf 

... 

y .  ^ 

♦  0 

5 

n  .  4  » / 1  4  ■.  n  ■*  \ 

-1 

1  .  2<* 

♦  1 

6 

-1  .  II. HI""  ■>><«  • 

*|. 

1  .  i/ft 

♦  1 

CfT|..n 

1 

3 . 0*.  2V/"W 

-1 

2.  2ft 

♦  1 

2 

3  .  o. 

-1 

2.  /ft 

♦  1 

1 

3  .f...  *  t> 

-1 

2.  2  ft 

♦  1 

2 

-3  .  Ai./bW"V  /* 

-1 

2.2ft 

♦  1 

i 

I.A.  7V  *>/ft 

-1 

2.2ft 

♦  1 

4 

3  .  ft.*’  'l  )•<•»«  2/ 

-1 

£  .  /ft 

•  1 

•> 

<*  SM.  •»*/*. 

-1 

1  ,i>4 

♦1 

6 

•  1  .  ."Oft 

•  U 

1  .Dl. 

♦  1 

ECT  I'IN 

1 

5 

-1  .  2S  t>Mi  .ftlft 

*  1 

H  ,  lift 

-1 

2 

-1  m* 

♦  1 

’/  ,  ft  .< 

-1 

1 

-1  ./S  %«<n 

•  1 

&  .Oft 

-1 

2 

1  .  2‘.  <•*.»  ,HU 

♦  1 

ft  .  Oft 

-1 

3 

-1  .  2S  1«M'  iO\ 

•  1 

ft  .lift 

-1 

4 

-1  .24  5"  U  w<l« 

♦  ‘I 

O  .'Ml 

-1 

b 

A.ATnl/  .■>!»» 

-1 

1  .  V. 

♦  1 

o 

•  1  .  (Mill'll  Ml  1)4 

♦  II 

1 .  lift 

♦  1 

f  r  t  i  i’n 

4 

1 

H.t'4S»>  '44ft 

♦  II 

1  . 2  * 

♦  0 

2 

/.ATI.  2.  lift 

♦  .1 

1  .  5ft 

♦  0 

y 

K.  111,4*'  *«4r. 

♦  d 

1.2ft 

♦  ') 

£ 

-H  .U/,S«>JS4ft 

♦  ll 

1.24 

♦  l) 

s 

A  .  *I4S'/>.>S4». 

♦  *r 

1  .2ft 

♦  0 

4 

ft  .  »UV  A  ' *v  4 n 

•  i. 

1.24 

♦  0 

b 

S  .  MS«W4.’VJ/, 

-1 

1 . 7  ft 

♦  1 

6 

-1  .  (Mll)«)U  M'U# 

♦  U 

1  .  Oft 

♦  1 

en  in* 

S 

i 

* .  M.S2u"4  3» 

♦  ii 

1  .14 

*0 

:  £ 

M  .  1  '.‘fcfc'MACift 

♦  li 

1.2ft 

♦  y 

i 

A  ,fl<.S<:'-'4  3ft 

♦  0 

1.1ft 

♦  y 

t 

-:«.*4S2 

♦V 

1.1ft 

♦0 

3 

n.«4S2i/"4)ft 

♦  II 

1.1.4 

♦y 

4 

*.H4S2‘"'4  -.*■ 

♦  l> 

1.1ft 

♦  0 

i 

'/£-. 

-1 

J.llft 

♦  i 

6 

-1  .  Ii.'hujm  •  ii u * 

♦  l* 

1  ."ft 

♦i 

ftTl-i* 

t . 

H  .H.'41  51  34ft 

♦  «i 

1.1ft 

♦n 

£ 

*  (l 

1  .  5ft 

♦u 

1 

d  .#.'>0  5  .  34  * 

♦  .l 

1.1ft 

♦  0 

2 

-A.H.'O  3  |  Snft 

♦  » 

1  .IK 

♦u 

3 

M  .  M  .*  4  ■>  3  1  <  4 

1  .I* 

♦  o 

4 

d.0,'03  1  3*  - 

♦  H 

1.1ft 

♦  0 

b 

1  .iMT  M  !  ft  y 

-3 

V  .  2ft 

♦  i 

6 

»1  .  IIOlH'ti  . ft 

♦  It 

1  . 1  It. 

♦  i 

3 

is 


i 


1 


Fig  6  concld 


'  7 


Fig  7 


TYPE  NUMBERS  FOR  FILTER  ORDER,  FILTER  TYPE,  CIRCUIT  TYPE  AND  LOW,  HIGH 
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TYPE  OVERALL  GAIN 
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TYPE  FREQUENCY  OF  PASS  BAND  EDGE 

200.0 
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TYPE  CHIN  AND  CMAX 
+■  1.0X-1 I  I .0&-7 
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1  .  0&3  I  .0X6 

BUTTERWORTH  FILTER 
NORMALISED  ROOTS 
SECTION  +1 

+  1 .8477590650X0  +1.00000000X0 

QUALITY  +5.41 1961001 5X- 1 

RES  FREQ.  1.00000000X0 

SECTION  +2 

+7.6536686473X-1  +1.00000000X0 

QUALITY  + 1 . 3065629649X0 

RES  FREQ.  1.00000000X0 

DENORMALISED  ROOTS 
SECTION  +1 

+2.3219625217X3  +1.5791367042X6 

QUALITY  +5.41 1961001 5&- 1 

RES  FREQ.  1.2566370614X3 
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TYPE  OP  AMP  INPUT  IMPEDANCE,  OUTPUT  IMPEDANCE  AND  OPEN  LOOP  GAIN 
«-  1.0&6  300.0  1.08.5 

TYPE  NUMBER  OF  EXTRA  FREQUENCIES 
<-  0 

IMPEDANCES_AND_GAIN_AT_FREQUENCY_Oj:OgOOOO&0_HZ 


INPUT  AND  OUTPUT  IMPEDANCES  OF  SECTION  1 
7.9577471 562&5  9.95455502108,-3 
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MAXIMUM  VOLTAGE  (AT  END  OF  SECTION  0)  IS  1.000048,0 
VOLTAGE  SPREAD  IS  1.00004&0 
TRANSMISSION  A-MATRIX  FOR  COMPLETE  CIRCUIT 
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7.06995782204-1 

1.799974629342 

(«  -3.011663542340  DB) 

FURTHER  ANALYSIS  OF  THIS  CIRCUIT  REQUIRED? 

«-  YES 

SENSITIVITY  ANALYSIS  REQUIRED? 
v  Y 

ARE  DIFFERENTIAL  SENSITIVITIES  REQUIRED? 

«-  N 

TYPE  NUMBER  OF  LARGE  CHANGES 
8 


Fig  7  contd 
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Fig  7  concld 


TYPE  TRIPLES  SPECIFYING  CHANGES 

Cl  2  2.0 

^  C2  2  2.0 

•*-  R1  2  5.0 

■*-  R2  2  5.0 

+  R3  2  5.0 

R4  2  5.0 

R5  2  5.0 

+  R6  2  2.0 

PERCENTAGE  CHANGE  IN  GAIN  AT  ZERO  FREQUENCY 
-1 .9607843148&0 

PERCENTAGE  CHANGE  IN  GAIN  AT  CUT-OFF  FREQUENCY 
-1 .9607843136&0 
IMPEDANCES  AND  GAINS  REQUIRED? 

«-  N 

FURTHER  ANALYSIS  OF  THIS  CIRCUIT  REQUIRED? 

•*-  H 

COMPONENTS  OR  RELATIVE  GAINS  TO  BE  ADJUSTED? 

+■  N 

IS  THERE  MORE  DATA? 

«-  N 


Fig  7  concld 


Fig  19 


a. 


4.  ■*:' 


Fig  22  Flowchart  for  data 


Fig  22  contd 
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